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1. Summary of Accomplishments 


The research objectives in this proposal were part of a continuing program at UW-CIMSS to 
develop and refine an automated geostationary satellite winds processing system which can 
be utilized in both research and operational environments. The majority of the originally 
proposed tasks were successfully accomplished, and in some cases the progress exceeded 
the original goals. Much of the research and development supported by this grant resulted in 
upgrades and modifications to the existing automated satellite winds tracking algorithm. These 
modifications were put to the test through case study demonstrations and numerical model 
impact studies. After being successfully demonstrated, the modifications and upgrades were 
implemented into the NESDIS algorithms in Washington DC, and have become part of the 
operational support. 

A major focus of the research supported under this grant attended to the continued 
development of water vapor tracked winds from geostationary observations. The fully 
automated UW-CIMSS tracking algorithm has been tuned to provide complete upper- 
tropospheric coverage from this data source, with data set quality close to that of operational 
cloud motion winds. Multispectral water vapor observations were collected and processed from 
several different geostationary satellites. The tracking and quality control algorithms were 
tuned and refined based on ground-truth comparisons and case studies involving impact on 
numerical model analyses and forecasts. The results have shown the water vapor motion 
winds are of good quality, compliment the cloud motion wind data, and can have a positive 
impact in NWP on many meteorological scales. 


2. Project Objectives and Results 


The original project objectives pertained to the improvement and upgrade of the existing 
UW-CIMSS automated satellite winds algorithm with a strong focus on attainment of high- 
quality wind information from the moisture channel observations. The results of these efforts 
are summarized here, with references given to publications containing further details. 


A) Improvements in tracer selection, height assignment methods, 
tracking and quality control strategies. 


The first step in obtaining high-quality wind vectors pertains to the identification of targets, 
or tracer selection. Traditional methods locate the highest pixel brightness values within a pre- 
defined target box and compute local gradients around those locations. Any of those gradients 
passing a pre-set threshold, and within a specified distance from a neighbor target, are 
selected as tracer targets. In reality, however, the optimum targets do not always correspond 
to the brightest pixel regions. An alternative strategy was explored, and eventually adopted. 
This new method reverses the above strategy. Two-dimensional gradients are computed for 
each pixel in the target box and the maximum gradients are identified. These gradients 
undergo a pixel brightness check to assure that the targets are clouds, and to filter clouds of 
certain levels (optional). This methodology insures the best targets (as defined by maximum 
gradients) are being retained in the tracer selection stage. 

Our targeting scheme also employs a spatial-coherence analysis which is applied to all 
prospective targets (Coakley and Bretherton, 1982). Brightness temperature data within the 
target domain is divided into as many 3X3 pixel boxes as possible, and means and standard 
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deviations of the 9 pixels within these boxes are computed. Boxes having standard deviations 
below a specified threshold are deemed to be “coherent” and are input into a 2-dimensional 
cluster analysis scheme to obtain the mean brightness for any coherent clusters (signals). Two 
checks are employed to filter out undesirable targets: 1) Only two coherent signals are allowed 
per target domain. The presence of more than 20 percent of coherent 3X3 boxes outside of 
the two largest clusters represents a multi-deck scene in most cases. Cloud-height diagnosis 
for multi-deck scenes is extremely difficult (Menzel et al., 1983), so these clusters are 
eliminated (area of ongoing research). 2) Scenes are not allowed to be too coherent. If more 
than 80 percent of the total number of 3X3 boxes within the target domain meet the standard 
deviation threshold, the scene is deemed to be too coherent and it fails as a potential target. 
These scenes typically represent uniform coastal features. 

As a result of the spatial-coherence filtering, the total number of targeting attempts within a 
given image can be increased significantly without any increase in the overall computational 
time necessary to process winds. The result is an efficient way to increase target density. The 
extra time used to complete a spatial-coherence analysis for all prospective targets is more 
than offset by the time saved by not having to height assign and track undesirable targets. 

The greatest contribution to vector RMS error pertaining to satellite-derived winds is 
attributed to assigning heights to the tracers. In regards to tracking features in the infrared, the 
development of the COrslicing algorithm (Menzel et al. 1983) was a major step forward in 
alleviating some of the height assignment problems associated with variable cloud emissivity 
(such as in cirrus). This is now the preferred method for assigning heights to semi-transparent 
clouds. However, the lack of a CO 2 absorption channel near 13 microns on most of the current 
geostationary satellites has necessitated the development of an alternative height assignment 
method. 
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The new method, termed the H 2 0-intercept algorithm (Nieman et al., 1993), is predicated 
on the fact that the radiances for two spectral bands vary linearly with cloud amount. 

Radiances from the infrared window and H 2 0 absorption band (6-7 microns) are measured 
and compared to calculations of radiances for both of these bands in opaque clouds at varying 
levels in an atmosphere specified by a numerical model prediction of temperature and humidity 
profiles. Measured and calculated radiances will agree for clear-sky and opaque cloud 
conditions. The cloud-top height is inferred from the linear extrapolation of measured 
radiances onto the opaque cloud calculations. The measured radiances used to infer the linear 
relationship between the two channels are the average radiances for the cluster of clearest 
(warmest) fields of view and the cluster of cloudiest (coldest) fields of view within the target 
area. If the calculated water-vapor radiances for clear sky is less than the measured water- 
vapor radiance, the calculated water-vapor radiances are adjusted to agree with the measured 
and the difference is attributed to an inaccurate transmittance used in the computation of the 
dear-column radiance. Calculated warm radiances that are greater than measured are not 
adjusted since the low measurement may be the result of doud contamination. 

Our evaluation of the H 2 0 intercept height assignment method has demonstrated that it is 
close in accuracy to the C0 2 slicing technique, and can be considered an acceptable 
replacement (Nieman et al., 1993). Our experience indicates that vector height assignments 
derived from this method yields standard deviations of 50-75 Mb when compared to 
rawinsonde ‘best fit’ values. This method is now being used in the UW-CIMSS algorithm with 
GOES-8 data, and has been incorporated into the algorithm used operationally by NESDIS. 

The second major step involved in satellite-derived wind vector production is target tracking. 
Standard correlation surface methods are employed in the UW-CIMSS algorithm. Image 
navigation and co-registration with subsequent images in the loop is critical to the accurate 
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determination of target displacement. As GOES-7 station keeping became less reliable, and 
the navigation procedures associated with the new three-axis stabilized GOES-8 were being 
fine tuned, it became apparent that the navigation strategy in the wind calculation code itself 
could be modified to improve the determination of target displacement and vector accuracy. 
Previously, it was assumed that all images were similarly navigated. The new software treats 
the navigation problem individually for each image of the loop using landmark techniques. In 
many cases, mis-registered images can be automatically corrected. This increases CPU time, 
but substantially reduces a source of vector error and improves the vector field coverage. 

The final important step in the process involves quality assurance. The UW-CIMSS 
algorithm incorporates an elaborate internal quality control (QC) system which is both objective 
and interactive. Certain tolerances and quality thresholds in the system are determined by the 
observing and processing strategy aimed at a particular application. Other controls are pre-set 
and based on empirically-determined thresholds. The QC system involves two major 
procedures. First, logic is entrained in the vector processing component of the winds algorithm 
and consists of a series of checks (Merrill et al. 1991). This step essentially weeds out 
unacceptable vectors based on factors such as low tracer correlation matches or high 
accelerations between vector pairs (3 image calculation). 

The second QC step involves an important and unique component of the UW-CIMSS winds 
processing algorithm, and is automatically invoked following the completed computation of the 
vector field. The basis of this procedure is a 3-dimensional recursive filter objective analysis 
with tunable interdependent parameters and tolerances, and an optional vector height re- 
assignment component (Hayden and Purser 1995). Numerous improvements have been 
incorporated into the latest version including: 1) tuning of the penalty function associated with 
the height re-assignment module to assign weights based on a minimization of vector RMS 
errors compared to rawinsonde ground truth, 2) addition of a maximum-allowable height re- 
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assignment, 3) incorporation of an empirically-based bias adjustment, which increments each 
vector with 7% of the forecast speed interpolated to the re-assigned height, 4) a tightened 
tropopause check to a threshold lapse rate of 0.5K per 25 Mb above 300 Mb to discourage 
stratospheric height assignment values, and 5) additional flexibility has been added to the 
code to allow for improved inspection and evaluation. 

B) Evaluation of cloud motion winds derived from the UW-CIMSS 
processing system. 

The aforementioned upgrades to the UW-CIMSS satellite winds processing system have 
resulted in improved vector fields in terms of coverage and accuracy. The assessment of these 
wind fields includes traditional comparisons with rawinsondes, and numerical model forecast 
impact studies. In terms of upper-level cloud motion winds produced by NESDIS (using the 
CIMSS algorithm) and determined from 30 minute GOES-7 IR imagery, verification statistics 
with collocated rawinsondes show a general improvement in vector bias and RMS error over a 
32 month period ending in July of 1995 (Nieman et al. 1996 and Appendix 6). Generally 
speaking, vector RMS error is around 7 m/sec, with a slow bias of less than 2 m/sec. Several 
of the upgrades mentioned above were implemented over the last 10 months of the period and 
resulted in noticeable improvements. For example, the bias adjustment in the objective editor 
became operational in December 1994 with a notable drop in the overall vector slow bias by 
almost 1 m/sec. 

Routine production of cloud motion winds by NESDIS (using the upgraded UW-CIMSS 
code) from GOES-8 IR imagery commenced in October of 1994. In general, preliminary 
verification statistics (Appendix 6) are promising and show the vector RMS error and bias to be 
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superior to the GOES-7 product by nearly .5 m/sec. Included in the processing algorithm for 
the GOES-8 winds are the latest tracer selection procedures, the H 2 0 intercept height 
assignment method in place of the C0 2 technique, and the upgraded objective editor including 
a 7% bias adjustment. 

More elaborate verification procedures involve the diagnosis of the impact of cloud motion 
wind data on numerical forecasts. As of this writing, two experiments have just been 
completed. The first involved NMC, and the impact of several weeks of GOES-8 cloud motion 
winds on their global assimilation system. The results were mixed. Impacts were found to be 
modest and in some cases negligible. Several factors may have led to this result, which is still 
being evaluated. The winds were processed at low density. Research (Velden et al. 1992) and 
recent experiments suggest higher density coverage may improve forecasts in general (more 
on this below). In addition, the impact was verified using 72-hour forecast global parameters. 
Since the winds were processed over a limited domain, these parameters may not reveal the 
true impact of the winds over the region which was covered. Further evaluation of this trial is 
anticipated. 

A second model impact experiment involved the processing of GOES-8 cloud motion winds 
at UW-CIMSS for 17 days in August at both high and low density (horizontal coverage). The 
two sets of winds were included into the UW-CIMSS Regional Assimilation System (CRAS). 
Forecast results were evaluated against rawinsondes and verifying analyses. In general, 
incorporation of the GOES-8 wind information at both densities resulted in improved CRAS 
forecasts of upper-level parameters (over the control runs which included only operationally 
available data only). Both skill scores and RMS error parameters were considered. Slightly 
improved impact results were realized with the higher-density data sets. These findings are 
being further evaluated and documented. 
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C) Development and evaluation of water vapor motion winds. 


Extraction of quantitative wind information from sequences of water vapor channel imagery 
is highly desirable to compliment information provided by cloud-tracked winds. A major focus of 
research under this contract was aimed at automating and optimizing the extraction of wind 
information from this data source. In this regard, significant advances were achieved. Our 
study has demonstrated that fully automated data sets which yield excellent horizontal 
coverage in the upper-troposphere are now achievable over a domain and on a time scale 
commensurate with real time operational constraints. In addition, results from real-time 
exercises show a high degree of data production success, indicating the reliability and 
robustness of the algorithm. These data sets can be produced on a medium-sized desktop 
computer. The UW-CIMSS algorithm has been tested using water vapor data from multiple 
geostationary platforms including GOES-7/8/9, Meteosat 3/5 and recently GMS-5. With the 
emergence of new-generation sensors (i.e., GOES-8/9, Meteosat-5) and optimized processing 
strategies, the quality of the product has improved to the point of being on par with operational 
upper-level cloud-tracked winds. Our processing strategy and major findings are summarized 
as follows: 

1) Water vapor winds - processing strategy 
a) A new target selection algorithm was developed that selects based upon the maximum 
gradient within the target box. Water vapor motion tracers in clear air are not primarily related 
to pixel brightness maxima (as in cloud tracking). Using gradient targeting, an increase in 
targets and resulting vectors of around 25% was found, with a .4 m/sec reduction in vector 
RMS error compared to collocated rawinsondes. 
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b) Image time separation is not as crucial compared to IR cloud motion wind tracking. 
However, 30 minute image intervals were found to be slightly superior to 1 hourly intervals, and 
are used routinely by UW-CIMSS in current satellite winds processing from the GOES imager 
data. 

c) Unlike cloud motion winds, uniform coverage is obtainable from water vapor vectors, 
and this folds into the quality control strategy. Generally speaking, certain constraints can be 
loosened (relative to the processing strategy with cloud motion winds) due to the plethora of 
‘neighbor checks’. Acceleration tolerances, deviations from the first guess, and the objective 
editor analysis ‘fit’ are good examples. Given the density and uniform coverage of vectors, 
questionable individual vectors are more easily identified and flagged as suspect or eliminated. 
On the other hand, fleets of ‘good’ vectors showing flow in disagreement with the first guess 
field have an improved chance of passing quality assurance criteria and impacting model 
analyses. 

2) Water vapor winds - characteristics and quality 

a) In terms of horizontal coverage in the upper-troposphere, water vapor vector fields are 
superior to cloud-motion derived fields due to the uniform nature of the coverage which allows 
delineation of meteorological features in clear air regimes (cf; Appendices 1, 2, 3). 

b) Mean vector field heights resulting from data obtained from different platforms will vary 
slightly depending on the attendant spectral response functions associated with each water 
vapor sensor. Generally speaking, vectors obtained from observations in the 6.5 - 7.0 micron 
region are assigned heights in the 150-500 Mb range, with a peak near 300 Mb. This is further 
illustrated and discussed in the 18-month Progress Report. 
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c) Using a large sample of GOES-8 vectors as an example, 47% of the original targets 
were traceable over time to yield wind vectors which passed correlation and acceleration 
constraints. The effects of the objective editor quality control were to reduce this yield by 
around 40%. Varying the weight given to the background guess at this stage had minimal 
effects on the overall quality. Down-weighting the guess by a factor of two resulted in 3% more 
vectors, many of which are found in active meteorological regimes such as circulations or high- 
amplitude troughs. 

d) Comparisons of vector quality (after objective editing) relative to collocated rawinsondes 
reveal a negligible speed bias, and vector RMS error values in the 6-8 m/sec range, which are 
comparable to operational upper-level cloud drift wind accuracies (cf; Appendices 1, 2, 3, 6). 
Relative to GOES-7, overall vector quality improves significantly (greater than .5 m/sec 
reduction in RMS error) when GOES-8 data is employed (Appendix 3). For model assimilation 
purposes, the variability of vector RMS error with wind speed (and indirectly, height) can 
provide important information. This relationship is illustrated in Fig. 1 , and shows the nearly 
linear increase in RMS error with vector wind speed. 

e) The quality of the data sets can also be measured in terms of their impact on numerical 
weather prediction. Several model impact studies conducted by UW-CIMSS and other 
modelling centers consistently show a modest positive impact on numerical predictions. 
Demonstration periods and case studies were used to show the impact of the winds on 
analyses and forecasts over scales ranging from global to sub-synoptic (hurricane). The 
specific results of these studies are outlined in Appendix 10. Modelling centers such as GFDL, 
NCEP, FSU and NRL are continuing to collaborate with us on the optimization of the water 
vapor winds assimilation towards a goal of maximizing numerical forecast improvement. 
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Figure 1. (Top) Accuracy of 6.7 micron GOES-8 water vapor motion vectors produced 
quasi— operationally by NESDIS compared to collocated rawinsondes (after objective 
QC), with respect to vector wind speed (3 m/s bins). (Bottom) Number of wind/raob 
matches per wind speed bin. 




3) Water vapor winds - sounder channels 

To supplement the upper-tropospheric wind vectors attainable from the 6 - 7 micron region, 
an investigation was conducted into the employment of water vapor channels borne on the 
GOES sounder. The GOES-8/9 sounders contain channels near 7.0 and 7.3 microns. The 
spectral response functions at these wavelengths carry radiometric sensitivities which peak 
lower in the troposphere relative to the 6.7 micron channel. Utilized in combination with the 6.7 
micron channel, the sounder water vapor channels can provide vertical information of the wind 
profile in the upper half of the troposphere (Fig. 2). While the sounder field of views are limited 
to 10 km (versus 8km for the imager), the signal to noise is improved over that of the imager. 
The result is a presentation of features in the moisture field which is comparable to that of the 
imager. 

The UW-CIMSS satellite winds tracking algorithm was applied using the 7.0 and 7.3 micron 
data. The targeting, tracking, processing and QC strategies were basically the same as with 
the 6.7 micron imager data (independent calibration tables were utilized). To get an idea of 
absolute accuracies, a large data set covering a 3 week period was produced over the US 
rawinsonde region. Comparisons with collocated rawinsonde winds yields accuracies 
comparable with, or slightly better than the 6.7 micron channel (6-7 m/sec). The accuracies 
with respect to wind speed are shown in Figs. 3 and 4, and as with the imager, the RMS errors 
increase with vector wind speed. An example of the typical coverage achievable from the UW- 
CIMSS automated system from each of the two sounder channels is shown in Fig. 5. Very 
good coherence in the vector fields is noted, as well as the general differences in the assigned 
heights. 

The value of the additional information from the sounder channels is especially obvious in 
cloudless regimes, where the vertical wind profile in the upper-troposphere at a particular 
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Fiaure 3. (Top) Accuracy of 7.0 micron GOES-8 sounder water vapor motion vectors 
compared to collocated rawinsondes (after objective QC), with respect to vector 
wind speed (3 m/s bins). (Bottom) Number of wind/raob matches per wind speed 
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Figure 4. (Top) Accuracy of 7.3 micron GOES-8 sounder water vapor motion 
vectors compared to collocated rawinsondes (after objective QC) , with respect 
to vector wind speed (3 m/s bins). (Bottom) Number of wind/raob matches per 
wind speed bin. 
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Figure 5. Example of wind coverage that can be achieved from tracking features 
in the 7.0 and 7.3 micron channels from the GOES-8 sounder. Yellow (blue) 
vectors are higher (lower) than 400 Mb. 






location can be gleaned. In fact, the multi-spectral approach frequently yields wind profiles 
which exhibit considerable shear (both directional and speed). The major limitation of the 
'surrogate imageri approach is the geographical coverage. Currently, the GOES-8 sounder is 
constrained to a domain which covers the continental US and adjacent waters. The real need 
for this data is over oceanic regions nearly void of mid-level wind information. Lobbying efforts 
are underway to free up the GOES-8 sounder during major marine events such as hurricanes. 
Starting January, 1996, the GOES-9 sounder will be pointed over the central/eastem North 
Pacific and will provide excellent coverage for this product to be tested and evaluated in terms 
of model forecast impact. 

4) Water vapor winds - specific applications and case studies 

Wind information from geostationary satellite water vapor bands can now be extracted on a 
nearly global scale. Water vapor vectors are produced operationally by EUMETSAT from the 
Meteosat-5. UW-CIMSS has demonstrated this capability from GOES 7 and 8, which will also 
be applicable from GOES-9. Preliminary testing at UW-CIMSS using the new GMS-5 water 
vapor channel indicates comparable information can be obtained. The uniform nature of the 
water vapor wind fields and quality comparable to operational cloud drift winds should lead to 
an operational product which nearly covers the full globe in the very near future. It is believed 
the continuous assimilation of this wind information will lead to reductions in global model 
forecast parameters. 

Preliminary studies indicate the impact the data can have on regional weather forecasts is 
encouraging (Appendices 2, 3, 10; 12-month Progress Report). The specific application of the 
data to hurricanes is a prime example, and has received considerable attention in this study. 
During the 1995 Atlantic hurricane season, sets of water vapor motion vectors were produced 
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for selected storms. In these cases, the processing strategy focused on high-density coverage 
in an attempt to capture subtle features in the hurricane circulation and environment. A good 
example is shown in Figs. 6 and 7 which depict data sets produced during tropical cyclone 
Pablo. In Fig. 6, tremendous outflow is apparent in the wind field concurrent with the storm 
intensification at that time. However, the vector field denotes a hostile upper-level environment 
upstream from the storm, which is about to impede on the hurricane outflow region. This 
hostile environment is indicated by the band of strong westerly winds (50 - 60 kts) located to 
the west of the storm. Strong vertical shear is well known to disrupt hurricane circulations and 
the intensification process, and in severe cases can decouple the low-level vortex from the 
convective support and effectively destroy the energetics of the storm. That is precisely what 
appears to have happened in this case. 24 hours later (Fig. 7) the storm has passed through 
the strong shear environment and is rapidly dissipating to a tropical depression. The strong 
band of westerlies that sheared Pablo was not indicated in any of the operational data 
available to the global analyses at this time, and intensity forecasts reflected this. Many other 
examples of features affecting storm motion and intensity were captured in the water vapor 
wind fields during the hurricane season, and selected cases are being investigated through 
model impact studies. Collaborations with GFDL, FSU, NRL-MRY and NCEP will yield model 
forecast impact results in the near future and these will be documented and reported on in the 
literature. 

The areal extent (global) and uniform nature of the coverage provided by water vapor 
tracked winds in the upper-troposphere can also lend itself to climatic applications. As an 
example, full disk water vapor wind sets were produced on a routine basis from Meteosat-3 
(stationed at 75W) during the Atlantic hurricane season of 1994. The data sets from July, 
August and September were combined into a 3-monthly mean wind field. In conjunction with 
this, upper-tropospheric humidity fields derived from the water vapor radiances were also 
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Figure 6. High-density GOES-8 water vapor motion vectors (6.7 micron channel) during Tropical Cyclone Pablo on October 7, 
1995. Vector colors refer to assigned heights (Mb) with legend in upper right hand corner. 





Figure 7. High-density GOES-8 water vapor motion vectors (6.7 micron channel) during Tropical Cyclone Pablo on October 8, 
1995. Vector colors refer to assigned heights (Mb) with legend in upper right hand corner. 





processed during the same period and averaged into a 3-month mean field. A divergence field 
was calculated from the mean wind field and compared to the mean humidity field. The result 
is shown in Fig. 8. A good match exists between dry (moist) regions and upper level 
convergence (divergence), which suggests regions of upper-tropospheric subsidence 
(convection). Of particular interest is the large region of mean subsidence over the eastern 
North Atlantic, which lies over the primary genesis area for Atlantic tropical cyclones. This data 
suggests the upper-level conditions over this region were not supportive for sustaining 
organized convection, which is the driving factor in hurricane genesis. In terms of tropical 
cyclone activity, the 1994 Atlantic hurricane season was well below normal. 

The fields depicted by the water vapor information may be providing a piece of evidence 
towards describing the inactive Atlantic hurricane season. The combination of the radiance 
data (transformed into upper-tropospheric humidity) and motion tracking (wind vectors) may 
also be useful in other climate diagnostic applications (Appendix 9). The implications are that 
this type of product may be an important addition to programs such as GEWEX or 
PATHFINDER. 
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Figure 8. Combination of 3-month mean upper-level tropospheric humidity and 
divergence derived from Meteosat-3 water vapor observations. The divergence 
field was computed from a 3-month mean water vapor motion vector field. 
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3. Future Directions 


The research and development supported under this grant is part of an ongoing effort to 
optimize the extraction of quantitative wind information from multispectral satellite 
observations. This general effort will continue under the ongoing support of the NOAA/NESDIS 
GOES evaluation program. The higher resolution and precision measurements afforded by the 
new generation of GOES satellites dictates the re-evaluation (upgrading) of sampling 
strategies. For example, it is likely we will be able to increase the density of vectors in the 
vicinity of meteorological events (i.e. storms, hurricanes) in order to capture the detail in the 
environmental wind field without sacrificing quality. Also, as processing capacity and 
assimilation systems advance, we foresee the processing and assimilation of satellite-derived 
data sets on an hourly basis in the near future. Further support for this R&D is being requested 
through the NASA Earth and Space Science (ESS) Project, in collaboration with NASA-GSFC 
scientists. 

The UW-CIMSS winds algorithm will also be applied and tuned to GMS-5 satellite 
observations in support of a grant to study the impact of satellite data on western North Pacific 
typhoon analysis and forecasting (funded through the Naval Research Laboratory). The goal 
of this effort is to demonstrate a fully automated satellite winds retrieval algorithm able to 
capture high-density wind vector fields in the environment of typhoons. The newly launched 
GMS-5 contains a water vapor channel, and wind observations from this data will be extracted 
and applied through NWP impact studies. 
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ABSTRACT 


Motions deduced in animated water vapor imagery from geostationary satellites can be used to infer wind fields 
in cloudless regimes. For the past several years, CIMSS has been exploring this potentially important source of 
global-scale wind information. Recently, METEOSAT-3 data has become routinely available to both the U.S. 
operational and research community. Compared with die current GOES satellite, the METEOSAT has a superior 
resolution (5km vs. 16km) in its water vapor channel. Preliminary work at CIMSS has demonstrated that wind sets 
derived from METEOSAT water vapor imagery can provide important upper-tropospheric wind information in data 
void areas, and can positively impact numerical model guidance in meteorological applications. Specifically, 
hurricane track forecasts can be improved. Currently, we are exploring methods to further improve the derivation 
and quality of the water vapor wind sets. 


1. INTRODUCTION 

Water vapor imagery from GOES satellites has been available for over a decade. These data are used 
extensively, mainly in a qualitative mode, by forecasters in the United States (Weldon and Holmes 1991). Some 
attempts have been made at quantifying the data by tracking features in time sequences of the imagery (Stewart et 
al. 1985; Hayden and Stewart 1987). For a variety of reasons, applications of this approach have produced marginal 
results (Velden 1990). Recently, METEOSAT-3 (M-3) was repositioned at 75W by the European Space Agency, in 
order to provide complete coverage of the Atlantic Ocean. Data from this satellite are being transmitted to the U.S. 
for operational use. Compared with the GOES satellite, the M-3 has; a superior resolution and signal-to-noise ratio in 
its water vapor channel, which translates into improved automated tracking capabilities. 

During the Atlantic hurricane seasons of 1992 and 1993, tracking algorithms were applied to the M-3 water 
vapor data in order to evaluate the coverage, accuracy and model impact of the derived vectors. Data sets were 
produced during several tropical cyclone cases, including Hurricane Andrew. In this paper, the M-3 water vapor 
wind sets are assessed, and their impact on a hurricane track forecast model is examined. 

2. STATISTICAL EVALUATION OF THE M-3 WATER VAPOR WINDS 

During the spring of 1992, M-3 water vapor wind sets were routinely produced on a daily basis (around 
1200UTC) at CIMSS. The domain of the wind sets covered the eastern U.S. and the western North Atlantic ocean 
basin (Fig. 1). The purpose of this exercise was to assess the horizontal and vertical coverage of the winds, and 
evaluate their accuracy. M-3 water vapor targets were selectively chosen to be constrained to approximate cloud- 
free areas in the imagery. This was accomplished by activating an empirically-determined brightness temperature 


threshold value which was not to be exceeded in the target selection step of the wind-tracking algorithm. An 
objective quality control procedure is invoked once the wind vector field has been derived (Hayden and Velden 
1991). 

The data set presented in Fig. 1 is typical of the data sets produced during the exercise. From a purely qualitative 
point of view, the horizontal coverage of the water vapor wind vectors shown in Fig. 1 is quite good in comparison 
to conventional and cloud-drift wind observations routinely available over the western North Atlantic basin. It was 
found that the vertical distribution of the assigned vector pressure-heights in cloud-free areas (based on brightness 
temperature to radiosonde temperature comparisons) was typically in the range of 200-500mb, with a peak near 
350m b. It was also clearly demonstrated for future considerations that the water vapor wind sets could be created on 
McIDAS (or VDUC at NMC) in a time scale commensurate with real time operations. 



Fig. 1. Typical M-3 water vapor wind set coverage. 


Comparisons between the M-3 water vapor wind vectors and collocated rawinsondes (within 1 .0 degree) over the 
period April - June 1992 were compiled and are presented in Table 1. The winds are also evaluated against the 
collocated first guess forecast (in this case the Aviation model 12h forecast). Both vector speed bias and RMS were 
computed and compared. The data sets are also stratified to reflect comparisons with selected rawinsondes in 
relatively remote areas. 


Table 1, Statistical evaluation of M-3 water vapor motion winds from the period April- June, 1992. 


1 ) Versus eastern US/westem 2) Versus Bermuda rawinsonde 3) Versus Guadeloupe rawinsonde 
Atlantic rawinsondes (N=981) (N=31) (N=50) 


H20 Winds NMC Forecast H20 Winds NMC Forecast 


H2Q Winds NMC Forecast 


Speed bias 
(m/s) 

-0.9 

-1.7 

-0.8 

-11 

-0.7 

-1.9 

Vector RMS 
(m/s) 

7.8 

6.6 

6.9 

7.1 

6.6 

6.8 


Overall, the M-3 vector speed bias (wind vector minus rawinsonde) is -0.9 m/sec, which is just under 1 m /sec 
better than the first guess forecast. The RMS error, however, is 1.2 m/sec higher than the first guess. The superior 
performance of the first guess RMS error can be explained by the fact that most of the comparisons were over the 
eastern United States, an area where the model has been properly initialized with abundant rawinsonde data. 


Examination of the relatively remote Bermuda and Guadeloupe rawin sonde comparisons, however, show the RMS 
error of the M-3 winds slightly below that of the first guess. The speed bias is also much improved over that of the 
first guess. 

More recently, M-3 wind sets have been routinely produced on a daily basis by NESDIS in Washington DC since 
May 1993. The algorithms and tracking methods were transferred from CIMSS in order to test and evaluate them in 
an operational environment The target selection and thresholds, and quality control procedures are nearly identical 
to those utilized at CIMSS. Validation against collocated rawinsondcs (N=1471) indicates a vector RMS of 8.2m/s, 
and a speed bias of -1.1 m/s, agreeing closely with the results in Table 1. 

3. IMPACT ON NUMERICAL HURRICANE TRACK FORECASTS 

Another way of quantitatively evaluating new data types is through model impact studies. In our investigation, a 
hurricane track forecast model (VICBAR, DeMaria et al. 1992) is used to test the sensitivity of the water vapor wind 
data on numerical hurricane track forecasts. VICBAR is a nested, spectral barotropic model that has been run in 
near-real time at the NOAA/Atlantic Oceanographic and Meteorological Laboratory Hurricane Research Division 
(AOML-HRD) for the past few years. The initial condition for the forecast model is a vertically averaged (mass 
weighted) deep-layer mean wind over the 850-200mb depth of the troposphere, with an added tropical cyclone 
bogus. The barotropic forecast model uses the shallow-water equations, with the forecast storm track determined 
from the location of the relative vorticity maximum on the innermost model mesh. This forecasting system has been 
used to evaluate impacts of other data types (Franklin and DeMaria 1992; Velden et al. 1992) 

Wind sets were produced daily at 1200UTC during several Atlantic tropical cyclones in 1992 and 1993. For the 
model impact evaluation, forecasts were included based on the following criteria: 1) the storm was within the 
domain of the generated wind set, and 2) the storm intensity was tropical storm strength or better (>35kts). There 
were 32 forecast cases, from 10 different storms (including Hurricane Andrew in 1992) that met these criteria. 

Table 2 shows the results of forecasts which were initialized with operationally available data plus the water 
vapor wind set produced at CIMSS, compared to control runs which were initialized on operationally available data 
only. Considering all of the cases together, the results show modest improvements to the forecasts with the inclusion 
of the water vapor wind data. The improvement in mean forecast error (MFi in Table 2) ranges from 2.29k at 24h to 
6.8% at 72h. After an adjustment is made for the serial correlation between forecast cases (Franklin and DeMaria 
1992), none of the improvements were found to be significant at the 95% confidence level. Another measure of 
forecast impact is the frequency of improved forecasts (FIF), which simply shows the percentage of forecasts which 
resulted in some improvement when the water vapor data were included. The FIF at 72h indicates that 71% (15 out 
of 21) of the forecasts were improved. 


Table 2. Impact of METEOSAT water vapor winds on VICBAR tropical cyclone track forecasts. The following 
verification statistics are valid for forecasts which included the water vapor winds: number of forecasts 
(N) , effective number of independent forecasts (N*), and the mean forecast error (MFE) relative to best 
track verification. Also given are comparisons with the control forecasts: the mean forecast improvement 
(MFI), expressed in both kilometers and as a percent relative to the control forecast error, the standard 
deviation of the improvements (SOI), the number of improved forecasts (IP), the frequency of improved 
forecasts (PIP) expressed as a percent, and whether the forecast improvements are statistically 
significant at the 95% confidence level (SIG) . 


Forecast 


Interval 

00 

N 

N* 

MFE 

(km) 

MFI 

(km) 

MFI 

% 

SDI 

(km) 

IF 

# 

FIF 

% 

SIG 

(y/n) 

24 

32 

26.6 

149.1 

3.3 

2.2 

23.8 

19 

60 

n 

48 

27 

22.4 

338.4 

3.9 

1.1 

48.2 

16 

60 

n 

72 

21 

17.4 

347.0 

23.6 

6.8 

91.6 

15 

71 

n 


Several of the selected tropical cyclone cases were relatively close to the U.S. mainland. It is reasonable to 
assume that in these cases, the VICBAR model was relatively well initialized by the nearby conventional 
(operational) data base, limiting the potential for the satellite data to have a positive impact. On the other hand, 
storms well out to sea should make better candidates for forecast improvement. To test this hypothesis, the sample 
was stratified to include only those cases east of 70W. The results of these 23 cases are shown in Table 3. While 
still not statistically significant at the 95% confidence limit, the MFI percentages and FIF are a notable 
improvement over the statistics presented in Table 2. 


Table 3. Same as Table 2, except only for cases east of 70W. 
Forecast 


Interval 

(h) 

N 

N* 

MFE 

(km) 

MFI 

(km) 

MFI 

% 

SDI 

(km) 

IF 

# 

FIF 

% 

SIG 

(y/n) 

24 

23 

19 

136.4 

9.0 

6.6 

20.2 

16 

70 

n 

48 

20 

17.8 

335.3 

10.2 

4.3 

60.1 

13 

65 

n 

72 

18 

16.2 

362.3 

30.4 

8.4 

98.6 

13 

72 

n 


Using the stratified sample, it is of interest to examine the distribution of forecast differences between the runs 
made with the water vapor winds and the control runs in order to fully appreciate the FIF results. Table 4 shows the 
distribution of the 72h forecast differences. Of the eighteen cases, only two forecasts that included the water vapor 
winds were notably degraded (greater than 100km) relative to the control forecasts. One of these poor forecasts was 
a Hurricane Bonnie case, which seemed to result at least in part because Bonnie's track did not follow the deep layer 
flow used in the model forecast. Rather, it followed a shallow-layer flow, in a direction quite different from the 
deep-layer flow. The water vapor winds provide measurements in the 200-500mb layer, thus affecting the upper part 
of the deep layer mean wind flow field. In this case, the control analysis approximated the shallow-layer flow more 
closely, and as a consequence, the VICBAR model control forecast without the winds was closer to the observed 
storm track. The reason for the other poor forecast has not yet been identified. 

On a more positive note. Table 4 shows that thirteen of the eighteen 72h forecasts were improvements over the 
control forecasts, with five of them being notable improvements (defined here as greater than 100km). Three of the 
notably improved forecasts occurred during Hurricane Andrew’s interaction with an upper-level cyclonic circulation. 
An example is shown in Fig. 2, from 1200 UT 19 August 1992. At this time the upper-level circulation was situated 
to the north of Andrew and quite evident in the water vapor imagery. The derived water vapor winds captured the 
circulation, as shown in Fig. 1. During the next 72h, Fig. 2 shows that Andrew’s track was only slightly affected by 
the upper-level low, shifting it north -northwest for a short time, before Andrew escaped its influence and turned 
more to the west However, the VICBAR control forecast from this time recurved Andrew to the north and 
eventually to the northeast in response to the influence of the upper low on the deep layer mean steering flow. 
Although the turn back towards the west was not predicted, the VICBAR forecast made with water vapor wind data 
responded with much less curvature, and an improved longer-range forecast. 


Table 4. Distribution of 72h forecast errors relative to control forecast errors for cases east of 70W. 

Forecast differences (km) 

* * * * ********** * * * * 

I 1 1 1 1 

200 100 0 100 200 



4. WORK IN PROGRESS 

To this point, the CIMSS water vapor winds have been produced from automated techniques designed for cloud 
drift wind extraction. This includes image preprocessing, target selection, the tracking algorithm, height assignment, 
and quality control. Clearly, all five of these basic procedures can be refined to better suit the characteristics of the 
water vapor data. Efforts are underway at CIMSS to examine these five areas to better tune them to water vapor 
tracking, and some examples are briefly discussed below. 

a) Image preprocessing 

NASA-MSFC has developed image filtering techniques to help reduce the noise in the water vapor imagery 
before the target selection process is invoked. These techniques will be investigated in collaboration with MSFC. 

b) Target selection 

The current target selection procedure involves identifying the coldest pixel in the target box with an 
accompanying brightness temperature gradient which satisfies certain thresholds. This philosophy was chosen to 
maximize the likelihood of selecting cloud edges as targets, since these features represent good tracers over time in 
conventional IR imageiy. Since we are attempting to provide vectors in cloud-free areas when using the water vapor 
data, this philosphy may not be optimal. A method will be investigated that simply searches for the maximum 
brightness temperature gradient over the target selection box. The resultant targets will be evaluated against 
routinely produced targets. 

c) Tracking algorithm 

The current automated tracking algorithm employed at CIMSS (Merrill 1989) uses a correlation method to 
identify potential matches of target radiance features in subsequent imagery. Currently, a sequence of three images, 
each separated by one hour, are used to derive two vectors. The two vectors are then averaged to achieve a final 
vector which is used to represent the tracer. Preliminary investigations indicate the averaging philosphy may not be 
appropriate with water vapor tracking since, compared to cloud tracking, correlations are lower and noise is higher. 
An algorithm will be investigated that selects the vector (first, second or average) that best satisfies accompanying 
quality control checks. In addition, the use of a four image sequence instead of three will be studied. The four image 
sequence will yield three vectors, which would provide a better evaluation of tracer consistency. 

d) Height assignment 

The method currently employed at CIMSS simply matches the target brightness temperature with a first guess 
temperature profile (assumes a one to one relationship) and extracts a pressure height A statistical evaluation of the 
actual relationship may yield superior vector height assignments (Hayden and Stewait 1987). 



c) Quality control 


A two-step quality control procedure is used in the CIMSS tracking system. The first step involves consistency 
checks between the first and second vector. The second step is an elaborate automated editor which analyzes the 
wind fields in conjunction with the first guess, and outputs confidence weights on each vector. The auto editor also 
attempts to find the optimum pressure height assignment for each vector, and therefore may reassign the original 
height to one which better fits the first guess profile. Efforts are underway to refine this system to take into account 
the water vapor vector characteristics. 

To optimize the extraction of quality motion vectors from METEOSAT water vapor imagery, CIMSS is 
collaborating with NASA-MSFC and EUMETSAT, and will also entertain ideas from the workshop. Model impact 
experiments are planned in 1994 in coordination with NMC, NASA-MSFC, ECMWF and AOML-HRD. 

5. SUMMARY 

The recent availability of METEOSAT-3 data over the western Atlantic Ocean has led to an effort to extract 
quantitative information from the water vapor channel. Data sets containing vectors derived from animated M-3 
water vapor imagery were produced during 1992 and 1993 using the CIMSS automated wind derivation algorithm. 
A statistical evaluation of the vectors reveals that the water vapor winds (relative to collocated rawinsondes) show a 
reduced speed bias compared to the collocated first guess forecast values. The vector RMS errors are larger (by 
about 1 m/sec) than the first guess over the eastern U.S., but become slightly lower than the first guess at remote 
locations (e.g., Bermuda and Guadeloupe). 

The water vapor wind sets were also demonstrated to have a slightly positive impact on barotropic numerical 
hurricane track forecasts (VICBAR model), although the results were not statistically significant at the 95% level. 
Except for a couple of examples, the preliminary results seem to suggest that the positive impact is maximized on 
cases well offshore and away from data-rich regions, as would be expected from intuitive reasoning. Most of the 
cases near the U.S. coast showed negligible or slightly negative impact. The most notable forecast improvements 
occurred during Hurricane Andrew's interaction with an upper-level circulation, which was well-captured by the M- 
3 water vapor winds. Overall, 72h track forecasts were improved by 7%, while nearly 3/4 of the VICBAR forecasts 
showed some improvement with the inclusion of the water vapor winds into the initial analysis. 

Current and future efforts are being directed towards refining the CIMSS wind extraction and quality control 
procedures to better account for the water vapor data and resultant motion vector characteristics. It is envisioned 
that these data can compliment cloud drift wind data to provide complete global upper-tropospheric wind coverage. 
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1. INTRODUCTION 

Water vapor imagery from geostationary satellites has been 
available for over a decade. These data are used extensively by 
operational analysts and forecasters, mainly in a qualitative 
mode (Weldon and Holmes 1991). In addition to qualitative 
applications, motions deduced in animated water vapor imagery 
can be used to infer wind fields in cloudless regimes, thereby 
augmenting the information provided by doud-drift wind vectors. 
Early attempts at quantifying the data by tracking features in 
water vapor imagery met with modest success (Stewart et al. 
1985; Hayden and Stewart 1987). More recently, automated 
techniques have been developed and refined, and have resulted 
in upper-level wind observations comparable in quality to current 
operational doud- tracked winds (Laurent 1993). 

In a recent study by Velden et al. (1993) it was demonstrated 
that wind sets derived from Meteosat-3 (M-3) water vapor 
imagery can provide important environmental wind information in 
data void areas surrounding tropical cydones, and can 
positively impact objective track forecasts. M-3 was repositioned 
to 75W by the European Space Agency in 1992 in order to 
provide complete coverage of the Atlantic Ocean. Data from this 
satellite are being transmitted to the U.S. for operational use. 
Compared with the current GOES-7 (G-7) satellite (positioned 
near 112W), the M-3 water vapor channel contains a superior 
horizontal resolution (5 km vs. 16 km). 

In this paper, we examine wind sets derived using automated 
procedures from both GOES-7 and Meteosat-3 full disk water 
vapor imagery in order to assess this data as a potentially 
important source of terge-sca/e wind information. As part of a 
product demonstration, wind sets were produced twice a day at 
Cl MSS during a six-week period in March and April (1994). 
These data sets are assessed in terms of geographic coverage, 
statistical accuracy, and meteorological impact through 
preliminary results of numerical model forecast studies. 
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2. DATA SET DESCRIPTION AND EVALUATION 

During March and April of 1994, M-3 and G-7 water vapor 
wind sets were routinely produced twice a day (around 00 and 
1 2 UTC) at CIMSS. The wind sets were processed from full-disk 
imagery using fully automated procedures on a rise-6000 
workstation. [See Merrill et al. 1991 for further details on the 
CIMSS automated winds processing system]. As a final step, an 
objective quality control procedure is invoked once the wind 
vector field has been derived (Hayden and Velden 1991). 

The field of winds presented in Fig. 1 is typical of the data 
sets produced during the exercise. From a purely qualitative 
point of view, the horizontal coverage of the water vapor wind 
vectors shown in Fig. 1 is quite good in comparison to upper- 
level conventional and doud-drift wind observations routinely 
available over oceanic areas. The very high-density coverage in 
the center of the figure is due to overlap of satellite coverage. It 
is noteworthy that, in general, there is a good agreement 
between satellites. It was found that the vertical distribution of 
the assigned vector pressure-heights (Hayden and Stewart 
1987; Hayden and Velden 1991) was typically in the range of 
200-500mb, with a peak near 300mb. It was also dearly 
demonstrated for future considerations that the water vapor 
wind sets could be created on a time scale commensurate with 
real time operations. 

Statistical comparisons between the water vapor wind 
vectors and collocated western hemisphere rawinsondes (within 
2.0 deg.) are shown in Table 1. Both vector speed bias and 
RMS were computed and compared. The results show values 
which are close to those of current operational doud-drift winds. 


Table 1 . Statistical evaluation of water vapor motion winds 
produced at CIMSS in March 1994, vs collocated western 
hemisphere rawinsondes (N - no. of matches). 


GOES-7 Meteosat-3 

N-2491 N«2791 

Speed bias (m/s) -.71 -.46 

Vector RMS (m/s) 7.93 7.89 


Fig. 1. Combined water vapor winds coverage from GOES-7 and Meteosat-3. 


3. MODEL IMPACT EXPERIMENTS - PRELIMINARY . 

FINDINGS 

Another way of quantitatively evaluating new data types is 
through model impact studies. The wind sets derived at CIMSS 
were made available in near real-time to the modelling groups at 
NMC, ECMWF and CIMSS (for the model impact segment of 
this exercise, only winds derived over marine areas were made 
available). Both NMC and CIMSS produced parallel forecasts in 
near real-time, while ECMWF collected the data sets for future 
evaluation. NMC conducted the excerise on a T62 version of the 
global spectral model. The CIMSS model covers a regional 
domain at 150 km horizontal resolution and 39 vertical levels. As 
of this writing, model impact results were just becoming 
available, and the findings presented here are preliminary. 
Further and more complete results will be presented at the 
conference. 

Overall, the impact of the water vapor data on upper- 
tropospheric forecasts in general can best be described as 
modestly positive, and case dependent The CIMSS model 


yields fairly consistent -5% improvements in 300 mb height 
RMS values in the 24-72 hr forecast range when evaluated 
against a verifying analysis over a North American domain. The 
impact can be much greater in certain synoptic situations 
(example given below). Very limited and preliminary results from 
the NMC forecast comparison is yielding very slight 
improvements to mean and RMS values of 200 mb forecast (72 
hr) wind fields compared with verifying analyses on a global 
domain. The slight magnitude of the improvements may be an 
inherent manifestation of the fact that the verification is being 
done over the full globe, whereas the wind sets are limited to 
roughly 1/3 global coverage. 

The impact of the data on NWP is perhaps best illustrated on 
a case by case basis. For example, a strong and well-defined 
circulation was evident in water vapor image loops on 24 March 
1994 off the western coast of the U.S. This particular storm 
system became deadly as it crossed the U.S. mainland. 
Torrential rains and resulting mudslides hit southern California 
as the storm entered the west coast Two days later, heavy 
rains and flooding was observed over many portions of the 








lower Mtesissippi and Ohio valleys. On 27 March, a series of 
tornadoes associated with the storm system Killed over 40 
people and resulted in extensive damage. Wind vectors derived 
on 24 March from animated water vapor sequences describe 
th is circulation rather well (Fig. 2A), and contributed Important 
information to the upper-level analysis of this event (Fig. 2B). 
Most of the vectors shown in Fig. 2A were assigned pressure- 
heights between 300 and 450 mb. It should be noted that there 
was little in the way of upper-level cloudiness associated with 
this particular circulation, making it a prime candidate to show 
the complimentary nature of the water vapor motion winds to 
cloud-drift winds (there were virtually no operational doud drift 
winds defining the upper-level drcutation). 




Fig. 2. A) Plot of water vapor motion winds describing the upper- 
level drculation off the west coast of the U.S. on 24 March 1994. 
B) Cl MSS model wind analysis (EXP) of this feature at 400mb. 


The impact of the water vapor data on the CIMSS model 
forecast of this event can be assessed by examining the 
differences between a control run (CON) which did not 
incorporate the water vapor wind data, and an experimental run 
(EXP) which induded this information. An example is illustrated 
in Fig. 3, which shows 300 mb height difference fields (EXP 
minus CON) from model forecasts initialized at 12 UT 24 March. 
Within the model package, a 3-D variational wind-mass 
adjustment routine is activated in areas lacking in mass 
information (i.e. oceanic regions). This routine essentially 
spreads the relative wealth of wind information onto the mass 
field through gradient balance adjustments. At the initial time, a 
perturbation in the height difference field (+20/- 10 m couplet) is 
observable off the west coast in the vidnity of the disturbance. 
This pattern is indicative of a sharper (and more potent) short 
wave trough in the EXP analysis compared to the control. This 
perturbation is maintained in the 24 hr forecast difference fields, 
and actually amplifies in the longer forecast intervals. By 48 hrs, 
the EXP forecast heights are greater then 30 m lower in 
advance of the transient upper-level trough associated with this 
event, and once again indicative of a stronger (and/or 
deepening) system relative to the control. These forecast 
difference fields are validated by comparing the forecast height 
fields from both runs to verifying analyses valid at the forecast 
time. To isolate the forecast impact of this particular system, the 
validation is only done over a region affected by the storm (i.e. 
the difference perturbations shown in Fig. 3). Using these 
criteria, the EXP height RMS errors relative to the CON are 12% 
lower at 24 hrs, 14% lower at 48 hrs, and 22% lower at 72 hrs. 

These results suggest that the addition of the water vapor 
wind information to the analysis of the disturbance off of the 
west coast positively impacted the upper-level forecast of this 
particular event through 72 hrs. Other forecasted parameters 
associated with this storm were influenced as well. An 
examination of the 72 hr precipitation prognoses indicates that 
both the EXP and the CON forecast a major rain event from 
northern Mississippi through South Carolina (not shown). The 
EXP forecast, however, shifted the precipitation maxium slightly 
to the west and increased it by 13 mm relative to the CON 
forecast. Overall, in an RMS sense, verification from rain gauge 
reports at 12 UT 27 March indicates the EXP forecast was 
superior. The 78 mm EXP forecast precipitation maximum was 
dose to a raingauge report of 72 mm located just 100 km to the 
west. 

4. SUMMARY 

It has been demonstrated that motions deduced utilizing 
water vapor imagery from geostationary satellites can be 
translated into useful upper-tropospheric wind information over 
data void areas. Preliminary results from model impact studies 
are indicating that this information can positively influence and 
impact NWP on regional scales, and possibly on global scales 
as well. Current and future efforts are being directed towards 
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Fig. 3. 300 mb height difference (m) fields (EXP minus CON) from Cl MSS regional model runs originating at 12 UT 24 March 1994. 


refining the Cl MSS wind extraction and quality control 
procedures to better account for the water vapor data and 
resultant motion vector characteristics. With the inclusion of a 
water-vapor imager on the next Japanese GMS satellite, it is 
envisioned that these data can compliment cloud drift wind 
observations to provide nearly complete global upper- 
tropospheric wind coverage. 
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1. INTRODUCTION 

Recent studies have demonstrated the 
feasibility of deriving wind information from 
animated sequences of geostationary water 
vapor imagery in order to supplement 
conventional observations (Laurent 1993; Velden 
1993; Holmlund 1993). Further work has 
examined the impact of these data on numerical 
weather prediction (NWP), and results suggest 
useful and complimentary upper-tropospheric 
wind information can be provided by these 
observations (Velden et al. 1993, Velden et al. 
1994). 

Development of algorithms at Cl MSS to 
extract quantitative wind information from 
geostationary satellite water vapor data has 
progressed to the stage of quasi-operational 
testing and demonstrations. Full-disk wind sets 
derived from two different satellites (GOES-7 and 
Meteosat-3) were routinely produced on a twice- 
daily basis in March-April 1994, and delivered to 
NMC and ECMWF for assimilation into research 
versions of^ their respective global assimilation 
systems. In addition, the wind sets were included 
in real-time assimilations using the CIMSS 
regional-scale model to assess finer-scale 
impacts on NWP. Preliminary impact results 
included here are encouraging, and the full 
assessment along with case study evaluations 
will be presented at the conference. 

With the recent successful launch of GOES- 
8, analysts and forecasters in the U.S. can 
expect satellite observations and products of 

* Corresponding author address: Christopher S. 
Velden, Univ. of Wisconsin-SSEC, 1225 W. 
Dayton St., Madison, Wisconsin 53706. 


unprecedented quality. One big improvement will 
be realized in the widely-utilized water vapor 
channel. The combination of higher resolution 
and reduced signal-to-noise will lead to an 
improved qualitative presentation for subjective 
analyses, but also should translate into superior 
quantitative products such as water vapor 
tracked winds. At the time of writing, GOES-8 
imagery was just becoming available. First results 
of tracking GOES-8 water vapor features are 
reported on in this paper, with further findings and 
applications to NWP to be presented at the 
conference. 

2. BACKGROUND 

Wind vectors derived from water vapor 
imagery are generated following the same 
general philosophy as with cloud drift wind 
tracking. Targets (structures in the imagery) are 
identified and tracked from successive images to 
derive motions (vectors), which are then assigned 
heights based on a best match between 
radiances (brightness temperatures) and a first 
guess temperature profile. In most cases, this 
height assignment will be representative of the 
layer of moisture being tracked, since the 
individual pixel radiances are achieved through 
contributions of energy over the depth of the 
attendant weighting function. Most vectors are 
assigned tropospheric pressure-heights in the 
200-500 mb range, with a peak near 300-350 mb. 
As a final step, an objective quality control 
procedure is invoked once the wind vector field 
has been derived (Hayden and Velden 1991). 
Procedures to derive wind sets have become fully 
automated (Merrill et al. 1991; Holmlund 1993), 
and it has been demonstrated these data sets 


can be created on a time scale commensurate 
with real time operations. 

Evaluation of the water vapor winds has 
centered on two approaches; statistical validation 
against rawinsondes. and model impacts. This 
paper focuses on the latter. Earlier studies have 
shown the water vapor tracked wind estimates to 
be comparable in quality to upper-tropospheric 
operational cloud tracked winds (Laurent 1993), 
with RMS values of 6.5-8 m/sec, and speed 
biases of less than 1 m/sec compared to 
collocated rawinsondes (Velden 1993; Velden et 
al. 1994). It was also demonstrated by Velden et 
al. 1994 that horizontal coverage is 
complimentary to wind information obtained 
through cloud tracking, and there is good 
consistency between observations derived from 
different satellites. 

A first attempt to assess the potential 
importance of the water vapor tracked wind 
information on NWP was reported in Velden et al. 
1993. In their study, a barotropic hurricane track 
forecast model was used to examine the impact 
of the winds on the analyses of Atlantic hurricane 
steering currents, and subsequent track 
forecasts. Modest improvements in track 
forecasts (relative to operational runs) were noted 
with the inclusion of the water vapor winds into 
the deep layer mean wind analyses which 
initialized the model. This was demonstrated by 
the fact that eleven out of fourteen 72 hr 
forecasts were improved, with a decrease in the 
mean forecast error of 8.2%. 

3. IMPACT ON REGIONAL AND GLOBAL 
ANALYSES AND FORECASTS 

Wind sets produced by CIMSS for a trial 
period in March- April 1994 (as described in 
Velden et al. 1994) were transmitted to NMC for 
model impact tests. NMC conducted the exercise 
using a T62 version of the global spectral model 
(similar to T126 AVN/MRF). Forecasts produced 
from analyses initialized with the water vapor 
winds are compared to operational forecasts in 
order to assess model impact. Shown in Table 1 
are comparisons of results from 72-h forecasts 
for a 15-day subset of the trial period. The 
forecast fields are verified against model 
analyses. Two latitude-dependent domains are 
considered in order to examine regional impacts. 


Table 1 . Impact of water vapor motion wind data 
on 72-h forecasts from the NMC global spectral 
model during a 15-day trial period in March, 
1994. 

60N-60S 

fsp 



opn 

exp i 

exp/opn 

200mb vector rms (m/s) 

12.00 

11.86 

9/6 

200mb acpsi 

.809 

.815 

9/6 

200mb acv 

.741 

.742 

8/7 

20N- 

20S 



200mb vector rms 

10.51 

10.22 

9/6 

200mb acpsi 

.654 

.681 

12/3 

200mb acv 

.541 

.550 

8/7 


opn = operational runs 

exp = experimental runs with water vapor winds 
fsp = frequency of superior predictions 
acpsi = anomaly correlation of stream function 
acv = anomaly correlation of v-component 


Considering first the near-global domain, the 
impacts are not dramatic, as might be expected 
since the wind coverage over the mid latitudes is 
limited. It should also be noted that the wind sets 
produced by CIMSS cover only a portion of the 
western hemisphere, while the verification 
statistics are global, further limiting the potential 
impact. Despite these caveats, the impacts are 
slightly positive in all parameters examined. 9 of 
the 15 forecasts were improved, with a slight 
reduction in the overall 200 mb wind vector rms. 

Greater impact is seen in the tropical band, a 
region of sparse conventional data, and abundant 
water vapor wind coverage. While still modest, a 
more substantial gain in forecast accuracy is 
noted, especially considering the rotational 
component of the wind (psi). For this particular 
parameter, 12 of the 15 forecasts were improved. 
These results suggest the winds are contributing 
to a better definition of circulations or wave-like 
features in the ITCZ region. 

The water vapor wind sets were also provided 
to ECMWF for their evaluation. While 


acknowledging the impressive coverage of the 
winds (Graeme Kelly, personal communication), 
ECMWF had not yet accomplished a rigorous 
model impact test at the time of this writing. 
However, plans are to evaluate the impact of the 
winds using a new 3D-VAR version of their global 
model, and these results will be presented at the 
conference. 

It is of interest to also examine the impact of 
the winds on regional scales. To accomplish this, 
forecast impact tests were conducted using the 
CIMSS regional assimilation system (CRAS), 
which is an adaptation of the Australian Bureau of 
Meteorology operational regional forecast model 
(Leslie et al. 1985). The model analysis covers a 
North American and adjacent waters domain at 
150 km horizontal resolution, with 39 vertical 
levels. The forecast model operates with 20 
vertical levels. Within the model package, a 3-D 
variational wind-mass adjustment routine is 
activated in areas lacking in mass information 
(i.e. oceanic regions). This routine essentially 
spreads the relative wealth of wind information 
onto the mass field through gradient balance 
adjustments. Experimental forecasts (EXP) 
including the water vapor winds are compared to 
control runs (CON) which contain only 
operationally available data. 

An example of data impact on the model 
analysis and forecast fields is illustrated in Figs, t 
and 2, respectively. Shown in Fig. t is a plot of 
the 275-325 mb water vapor winds (large vectors) 
for 00 UTC on 22 March 1994. It should be noted 
that winds were only produced over marine 
areas. The vectors are overlain on a 300 mb 
model wind analysis difference field (EXP - 
CON), indicated by the smaller barbs (a barb 
length of one grid point interval equals 8 m/sec). 
As can be seen, differences at this particular 
level are modest, but can be as great as 6 m/sec 
such as in the trough off of the east coast of the 
US. 

These analysis differences translate into 48-h 
forecast differences shown in Fig. 2. The 48-h 
EXP forecast 300 mb height field is overlain to 
show the large-scale pattern. In order to highlight 
the differences, in this figure a barb length of one 
grid point interval equals 5 m/sec. Fig. 2 shows 
small differences cover most of the domain, with 
notable differences in several locations including 


the EXP forecast tendency to amplify the ridge 
over the central U S. 

Model impact results from CRAS forecasts of 
19 cases from the March/April 1994 trial period 
reveal a trend similar to that found with the NMC 
impact study. In general, a modest but positive 
improvement in RMS and bias statistics of upper- 
level height and wind fields are noted. An 
example of the forecast impact (72 hr) in terms of 
RMS at 400 mb is given in Fig. 3. In this figure, a 
comparison of CON and EXP forecast RMS 
errors (versus rawinsondes) is presented as a 
CON minus EXP value (meters) for each case. 
Positive model forecast impacts (reduced RMS 
errors as indicated by cases above the 0 line) 
produced by the addition of the water vapor 
winds to the analyses are indicated in 12 of the 
19 cases. Only 3 cases exhibit a degraded 
forecast, and 4 had virtually no impact. These 
results are fairly representative for mandatory 
levels between 200 and 500 mb, where the 
greatest impact would be expected due to the 
vertical distribution of the data. 

The impact of the winds on the regional 
model forecasts appears to be highly case 
dependent. While the absolute forecast error 
improvements are not large in a domain-wide 
RMS sense, regional forecasts may be more 
significantly affected. Qualitative examination of 
the cases run with the CRAS model indicates 
greatest positive impact occurs with active 
patterns in the eastern Pacific (i.e. upper-level 
circulations or jet streaks). These active periods 
correspond to some of the positive impact peaks 
in Fig. 3. Often these systems are not adequately 
defined by operationally available data (i.e. 
aircraft reports or cloud tracked winds). For an 
example, see Velden et al. (1994). 

4. PRELIMINARY DATA FROM GOES-8 

At the time of writing, data from GOES-8 were 
becoming available as part of system check-out. 
Although limited, observations were collected for 
several time periods in August 1994. 
Qualitatively, the GOES-8 water vapor imagery 
appears to be superior to the GOES-7 imagery, 
due to improved resolution (8 km vs. 16 km) and 
signal-to-noise. Wind algorithms were applied to 
extract water vapor vectors to assess (in a 
preliminary sense) and compare to GOES-7 



Fig. 1 . Plot of 275-325 mb water vapor winds (large barbs, full stick = 5 m/sec and flag = 25 m/sec) 
over 300 mb CRAS wind analysis difference (EXP minus CON) field (arrow length of one grid length 
interval = 8 m/sec) for 00 UTC 22 March 1994. 



Fig. 2. Difference field (EXP - CON) from two CRAS 48 hr forecasts of 300 mb wind (arrow length 
one grid interval * 5 m/sec) initialized at 00 UTC 22 March 1994. Also shown (contours in meters) 


the 300 mb EXP height forecast. 
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Fig. 3. Comparison of RMS errors (verified against rawinsondes) for CRAS model 72 hr forecasts of 
400 mb heights for 19 cases in March/April 1994, expressed in terms of RMS(CON) minus 
RMS(EXP). Positive impact of the water vapor winds (EXP) is indicated in cases with points above the 
zero line. A few of the days contain multiple data sets (00 and 12 UTC). 



Fig. 4. Example of the horizontal coverage of water vapor winds attainable from GOES-8, plotted in 
knots over a GOES-8 water vapor image. 



winds. An example of the typical coverage 
achieved with GOES-8 is shown in Fig. 4. Wind 
vector coverage is quite remarkable, and 
indicates a high degree of success in obtaining 
traceable features. 

A limited inter-comparison of vector quality 
(GOES-8 vs. GOES-7) is shown in Table 2. For 7 
days in August, collocated (space and time) wind 
vectors from each satellite were derived, quality- 
controlled and are compared to collocated 
rawinsondes. The results show a -10% reduction 
in GOES-8 wind vector RMSE relative to GOES- 
7. While preliminary and limited, these results are 
quite encouraging. This assessment will continue 
as GOES-8 data become routinely available, and 
an update will be presented at the conference. 


Table 2. Inter-comparison of GOES-8 and 
GOES-7 water vapor motion wind vectors verified 
against collocated rawinsondes for a 7 day period 
in August, 1994. NUM = number of comparisons. 
RMSE and BIAS are in m/sec. 



RMSE 

BIAS 

NUM 

GOES-7 

6.96 

0.10 

494 

GOES-8 

6.32 

0.09 

494 


5. SUMMARY 

Preliminary results from NWP impact studies 
are indicating that upper-tropospheric wind 
information provided by tracking motions in 
sequences of geostationary satellite water vapor 
imagery can positively influence forecasts on 
regional scales, and possibly on global scales as 
well. The data are complimentary to cloud- 
tracked winds by providing data in cloud-free 
regions, as well as comparable in quality. First 
results from GOES-8 winds are encouraging, and 
further efforts and model impacts will be directed 
towards optimizing these data in NWP. Assuming 
successful launches of GOES-J and GMS-5 
satellites in 1995, high quality and resolution 
water vapor imagers will be available to provide 
nearly complete global upper-tropospheric wind 
coverage. 
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1. INTRODUCTION 

Water vapor imagery from GOES satellites has been 
available for over a decade. These data are used 
extensively, mainly in a qualitative mode, by forecasters in 
the United States (Weldon and Holmes 1991). Some 
attempts have been made at quantifying the data by 
tracking features in time sequences of the imagery (Stewart 
et al. 1985; Hayden and Stewart 1987). For a variety of 
reasons, applications of this approach have produced 
marginal results (Velden 1990). Recently, METEOSAT-3 
(M-3) was repositioned at 50W by the European Space 
Agency, in order to provide complete coverage of the 
Atlantic Ocean. Data from this satellite are being 
transmitted to the U.S. for operational use. Compared with 
the GOES satellite, the M-3 has a superior resolution and 
signal-to-noise ratio in its water vapor channel, which 
translates into improved automated tracking capabilities. 

During a period in 1992 which included the Atlantic 
hurricane season, water vapor tracking algorithms were 
applied to the M-3 data in order to evaluate the coverage, 
accuracy and model impact of the derived vectors. Data 
sets were produced during several tropical cyclone cases, 
including Hurricane Andrew. In this paper, the M-3 water 
vapor wind sets are assessed, and their impact on a 
hurricane track forecast model is examined. 

2. STATISTICAL EVALUATION OF THE M-3 WINDS 

During the spring of 1992, M-3 wind sets were 
routinely produced on a daily basis (around 12 UTC) at 
CIMSS. The domain of the wind sets covered the eastern 
U.S. and the western North Atlantic ocean basin (Fig. 1). 
The purpose of this exercise was to assess the horizontal 
and vertical coverage of the winds, and evaluate their 
accuracy. M-3 water vapor targets were selectively chosen 
to be constrained to cloud-free areas in the imagery. This 
was accomplished by activating an empiricaDy-deternrined 


brightness temperature threshold value which was not to 
be exceeded in the target selection step of the wind- 
tracking algorithm. 

The data set presented in Fig. 1 is typical of the data 
sets produced during the exercise. From a purely 
qualitative point of view, the horizontal coverage of the 
wind vectors shown in Fig. 1 is quite good in comparison 
to conventional observations routinely available over the 
western North Atlantic basin. It was found that the vertical 
distribution of the assigned vector pressure-heights was 
typically in the range of 200-500mb, with a maximum near 
350mb. It was also clearly demonstrated for future 
considerations that the wind sets could be created on 
McIDAS (or VDUC at NMC) in a time scale 
commensurate with real time operations. 



Comparisons between the M-3 wind vectors and 
collocated rawinsondes (within 1 .0 degree) were compiled 
and are presented in Table 1 . The winds are also evaluated 
against the collocated first guess forecast (in this case the 
Aviation model 12h forecast). Both vector speed bias and 
RMS were computed and compared. The data sets are also 
stratified to reflect comparisons with selected rawinsondes 
in relatively remote areas. 


Table 1. Statistical evaluation of M-3 
water vapor motion winds 

1) Versus eastern US/western Atlantic rawinsondes 
(N-981) 

H20 Winds NMC Forecast 

Speed Bias (m/s) -0.9 -1.8 

Vector RMS (m/s) 6.5 5.5 

2) Versus Bermuda rawinsonde (N*31) 

H20 Winds NMC Forecast 

Speed Bias (m/s) -0.8 -2.1 

Vector RMS (m/s) 6.2 6.4 

3) Versus Guadeloupe rawinsonde (N=»50) 

H20 Winds NMC Forecast 

Speed Bias (m/s) -0.7 -1.9 

Vector RMS (m/s) 5.6 5.9 


Overall, the M-3 vector speed bias (wind vector minus 
rawinsonde) is -0.9 m/sec, which is about 1 m/sec better 
than the first guess forecast. The RMS error, however, is 
1 m/sec higher than the first guess. The superior 
performance of the first guess RMS error can be explained 
by the fact that most of the comparisons were over the 
eastern United States, an area where the model has been 
properly initialized with abundant rawinsonde data. 
Examination of the relatively remote Bermuda and 
Guadeloupe rawinsonde comparisons, however, show the 
RMS error of the M-3 winds slightly below that of the first 
guess. The speed bias is also much improved over that of 
the first guess. 

3. IMPACT ON NUMERICAL HURRICANE TRACK 
FORECASTS 

Another way of quantitatively evaluating new data 
types is through model impact studies. In our investigation, 
a hurricane track forecast model (VICBAR, DeMaria et al. 
1992) is used to test the sensitivity of the water vapor wind 
data on numerical hurricane track forecasts. VICBAR is a 
nested, spectral barotropic model that has been run in 
near-real time at the NO AA/ Atlantic Oceanographic and 
Meteorological Laboratory Hurricane Research Division 


(AOML-HRD) for the past few years. The initial condition 
for the forecast model is a vertically averaged (mass 
weighted) deep-layer mean wind over the 850-200mb 
depth of the troposphere, with an added tropical cyclone 
bogus. The barotropic forecast model uses the shallow- 
water equations, with the forecast storm track determined 
from the location of the relative vortici ty maximum on the 
innermost model mesh. This forecasting system has been 
used to evaluate impacts of other data types (Franklin and 
DeMaria 1992; Velden et al. 1992) 

Wind sets were produced daily at 12UTC during 
several Atlantic tropical cyclones in 1992. For the model 
impact evaluation, forecasts were included based on the 
following criteria: 1) storm within the domain of generated 
wind set, and 2) storm intensity of tropical storm strength 
or better. There were 19 forecast cases, from 4 different 
storms (including Hurricane Andrew) that met these 
criteria (Fig. 2). 



forecasts were made (circles) in this study. 

Table 2 shows the results of forecasts which were 
initialized with operationally available data plus the water 
vapor wind data, compared to control runs which woe 
initialized on operationally available data only. 
Considering all of the cases together, the results show 
modest improvements to the forecasts with the inclusion of 
the water vapor wind data. The improvement in mean 
forecast error (MFI in Table 2) ranges from 1.8% at 24h 
to 8.2% at 72h, however, a fter an adjustment is made for 
the serial correlation between forecast cases (Franklin and 
DeMaria 1992), none of the improvements were found to 
be significant at the 95% confidence level. Another 
measure of forecast impact is the frequency of improved 
forecasts (FIF), which simply shows the percentage of 
forecasts which resulted in some improvement when the 
water vapor data were included. The FIF at 72h indicates 
that 79% (1 1 out of 14) of the forecasts were improved. 




Table 2. Impact o f METEOSAT water vapor winds on 
VTCBAR tropical cyclone track forecasts. The 
following verification statistics are valid for 
forecasts which included the water vapor winds: 
number of forecasts (N), effective number of 
independent forecasts (N*), and the mean forecast 
error (KFE) relative to best track verification. 
Also given are comparisons with the control 
forecasts: the mean forecast improvement (MFI), 
expressed in both kilometers and as a percent 
relative to the control forecast error, the 
standard deviation of the improvements (SDI ) , the 
number of improved forecasts (IF), the frequency 
of improved forecasts (FIF) expressed as a 
percent, and whether the forecast improvements are 
statistically significant at the 95Z confidence 
level (SIG). 


Forecast 
Interval 
th) H 

N* 

MFE 

(km) 

MFI 

(km) 

MFI 

Z 

SDI 

(km) 

IF 

# 

FIF 

Z 

SIG 

(y/n) 

24 

19 

16 

139.9 

2.6 

1.8 

22.9 

10 

53 

n 

48 

17 

14.4 

354.1 

6.8 

1.9 

49.1 

11 

65 

n 

72 

14 

11.6 

381.2 

34.2 8.2 

94.5 

11 

79 

n 


From an examination of Fig. 2, it is evident that 
several of the selected cases are relatively close to the U.S. 
mainland. It is reasonable to assume that in these cases, 
the V1CBAR model was relatively well initialized by the 
nearby conventional (operational) data base, limiting the 
potential for the satellite data to have a positive impact. On 
the other hand, storms well out to sea should make better 
candidates for forecast improvement. To test this 
hypothesis, the sample was stratified to include only those 
cases east of 70W. The results of these 13 cases are shown 
in Table 3. While still not significant, the MFI per cen t a ges 


Table 3. Same as Table 2, except only cases east 


of 70W. 
Forecast 
Interval 
(h) H 

H* 

MFE 

(km) 

MFI 

(km) 

MFI 

Z 

SDI 

(km) 

IF FIF 

# z 

SIG 

(y/n) 

24 13 

10.8 

104.5 

8.9 

7.8 

18.8 

9 

69 

n 

48 13 

10.8 

221.6 

10.7 

4.6 

53.2 

10 

77 

n 

72 13 

10.8 

389.0 

35.9 

8.5 

107.2 

10 

77 

n 


are a notable improvement at 24 and 4£h over the sample 
presented in Table 2. The FIF are also improved. The 
values at 72h are nearly unchanged since only one case 
was deleted from the original sample. 

It is of interest to examine the distribution of forecast 
differences in the stratified sample, in order to fully 
appreciate the FIF results. Table 4 shows the distribution 
of the 72h forecast differences. Of the 13 cases, two 
forecasts were notably degraded by the inclusion of the 
water vapor winds. One of these poor forecasts was a 
Hurricane Bonnie case, which seemed to result at least in 
part because Bonnie's track did not follow the deep layer 
flow used in the model forecast. Rather, it followed a 
shallow-layer flow, in a direction quite different from the 
deep-layer flow. The water vapor winds provide 
measurements in the 200-500mb layer, thus affecting the 
upper part of the deep layer mean wind flow field. In this 
case, the control analysis approximated the shallow-layer 
flow more closely, and as a consequence, the VICBAR 
model control forecast without the winds was closer to the 
observed st or m track. The reason for the other poor 
forecast has not yet been identified. 


Table 4. Distribution of 72~h forecasts relative 
to control forecasts for cases east of 70W. 

Forecast differences (km) 

* * ** * *** * * * 


200 100 0 100 200 
BETTER MORSE 


On a more positive note, Table 4 shows that 10 of the 
13 72h forecasts were improvements over the control 
forecasts, with 4 of them being notable improvements 
(defined here as greater than 100km). Three of the notably 
improved forecasts occurred during Hurricane Andrew's 
interaction with an upper-level cyclonic circulation. An 
example is shown in Fig. 3, from 12 UT 19 August 1992. 
At this time the upper-level circulation was situated to the 
north of Andrew and quite evident in the water vapor 
imagery. The derived water vapor winds captured the 
circulation, as shown in Fig. 1. During the next 72h, Fig. 
3 shows that Andrew's track was only slightly affected by 
the upper-level low, shifting it NNW for a short time 
before Andrew escaped its influence and turned more to 
the west The VICBAR control forecast from this time 
recurved Andrew to the north and eventually to the 
northeast in response to the influence of the upper low on 
the deep layer mean steering flow. Although the turn to the 
west was not predicted, the VICBAR forecast made with 
water vapor wind data responded with much less 
curvature, and an improved longer-range forecast 
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Fig. 3. Example showing impact of water vapor 
winds on a VICBAR forecast of Andrew's track 


4. SUMMARY 

The recent availability of METEOSAT data over the 
western Atlantic Ocean has led to an effort to extract 
quantitative information from the water vapor channel. 
Data sets containing vectors derived from animated water 
vapor imagery were produced during 1992 using the 
CIMSS automated wind derivation algorithm. A statistical 
evaluation of the vectors reveals that the water vapor 
winds (relative to collocated rawinsondes) show a reduced 
speed bias compared to the collocated first guess forecast 
values. The vector RMS errors are larger (by about 1 
m/sec) than the first guess over the eastern U.S., but 
become slightly lower than the first guess at remote 
locations (e.g., Bermuda and Guadeloupe). 

The wind sets were also demonstrated to have a 
slightly positive impact on barotropic numerical hurricane 
track forecasts (VICBAR model), although the results were 
not statistically significant at the 95% level. Except for a 
couple of examples, the preliminary results seem to 
suggest that the positive impact is maximized on cases well 
offshore and away from data-rich regions, as would be 
expected from intuitive reasoning. Most of the cases near 
the U.S. coast showed negligible or slightly negative 
impact. The most notable forecast improvements occurred 
during Hurricane Andrew's interaction with an upper-level 
low, which was well-captured by the water vapor winds. 
Overall, 72h track forecasts were improved by an average 
of around 8%, while nearly 80% of the VICBAR forecasts 
showed some improvement with the inclusion of the water 
vapor winds into the initial analysis. 
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3.5 APPLICATIONS OF GOES-8/9 DATA TO HURRICANE ANALYSIS 


Christopher S. Velden 

Cooperative Institute for Meteorological Satellite Studies 
University of Wisconsin, Madison, Wisconsin 


1. INTRODUCTION 

GOES-8 was successfully launched in 1994, 
and has performed up to expectations and in 
some cases has actually exceeded pre-launch 
specifications. A similar geostationary satellite, 
GOES-9, was launched in May of 1995. This 
new generation of NOAA’s geostationary 
satellites carry a superior design and 
instrumentation package that allow for greater 
detection of meteorological features and 
parameters. The new GOES imager has a 5 
band multispectral capability with high spatial 
resolution, while the sounder contains 18 
thermal infrared (IR) bands plus a low- 
resolution visible band. The imager carries a 
visible channel with 10-bit quantization and 
increased sampling frequency, a short-wave 
and long-wave window channel, and a water 
vapor band with a twofold increase in spatial 
resolution and a factor of 3 improvement in 
signal-to-noise over that obtained from previous 
GOES sensors. 

The advances in observing the earth's 
atmospheric system anticipated from these 
improvements are outlined in Menzel and 
Purdom (1994). The specific impact of this 
improved remote sensing capability on the 
analysis of tropical cyclones is discussed here. 

2. HURRICANE ANALYSIS 

The detector improvements summarized 
above have lead to superior qualitative image 
presentations of hurricane structure and 
development. In the visible spectrum, sharper 
images are resulting in improved cloud top and 
cloud edge feature detection. This helps in the 
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subjective interpretation of the imagery in 
regards to eye and eyewall characteristics, but 
also enhances the ability of automated cloud 
tracking algorithms. In addition, with proper 
image display and enhancement techniques, 
the 10-bit imagery allows extended use in low 
light situations to aid in eye location. 

In the IR spectrum, the improved spatial 
resolution and signal-to-noise again benefits the 
ability to detect trackable cloud features, and is 
resulting in higher quantity and quality wind 
vectors. Velden et al. (1992) have shown that 
high-density satellite-derived wind information 
in the hurricane near-environment can lead to 
important forecast improvements. Rainfall 
estimation techniques should also benefit from 
the higher quality data. Determination of sea 
surface temperatures (important to hurricane 
energetics and potential intensity) by split 
window techniques should become more 
accurate, especially with the availability of the 
short wave window channel. 

Water vapor observations from geostationary 
satellites have been shown to be useful in 
delineating features in the hurricane 
environment that can affect intensity and 
motion (Dvorak 1984; Weldon and Holmes 
1991, Velden et al., 1993). Upper-level wind 
data sets derived utilizing water vapor image 
sequences from previous geostationary 
satellites have resulted in modest positive 
impact on numerical weather prediction (Velden 
1995a). The water vapor images from GOES-8 
are of unprecedented quality, and indications 
(Velden 1995b) are that image sequences are 
providing substantially improved wind vectors in 
cloud-free regimes (more on this below). 

The importance of the sounder information 
is to provide multivariate vertical profile 
observations in doud-free regions in the 
hurricane environment, especially valuable in 
regions void of conventional data (i.e. 
subtropical ridge in the North Atlantic). Thermal 
gradient winds derived from these profiles have 
been shown to benefit hurricane track forecasts 


(Velden et al. 1992). Preliminary indications are 
that the improvements in the GOES-8/9 
sounder detection capabilities are resulting in 
superior gradient wind information. 

As mentioned above, quantitative wind 
information derived from sequences of water 
vapor imagery has the potential to impact 
hurricane analysis by improving the depiction of 
the upper-level environmental flow fields. 
During hurricanes Florence and Gordon (1994), 
data sets were produced at 12-hour intervals 
and delivered to numerical modeling centers for 
their evaluation on track forecast impacts. First 
results obtained from the Hurricane Research 
Division’s barotropic hurricane forecast model 
(V1CBAR) were encouraging. The wind sets 
were shown to yield positive impact on 72-hour 
track forecasts (nearly 2/3 of the forecasts were 
improvements over the control runs). However, 
the improvements were quite modest, which 
might be partially explained by the 
characteristics of the VICBAR model, which 
initializes on a deep layer mean wind analysis 
that gives little weight to the upper-levels, and 
contains no vertical blending of data above 
350mb (most of the water vapor vectors are 
assigned heights above 350mb). 

For this reason, the data sets are being 
tested in more sophisticated hurricane track 
assimilation and forecast systems (FSU, GFDL) 
which employ physical initialization schemes 
and use the full primitive equations. The results 
of these trials will be presented at the 
conference. 


3. SUMMARY 

GOES-8 has performed up to expectations, 
and in some cases is exceeding pre-launch 
specifications. Applications of this new and 
improved satellite data to the analysis of 
hurricanes commenced in 1994 during the 
spacecraft and instrumentation check-out, and 
continued during the active 1995 Atlantic 
hurricane season. Findings are revealing that 
the data being provided by GOES-8 are 
superior to observations from previous 
geostationary satellites, and have the potential 
to positively impact both the subjective and 
objective analysis of hurricane situations. 
Development and testing of algorithms to 
optimally extract useful information from the 
observations is underway. 

GOES-9 was launched in May of 1995, and 
positioned at 90W during the 4-6 month 


checkout phase. This allowed dual-satellite 
coverage and inter-satellite data comparisons. 
Preliminary indications are that the GOES-9 
observations are as good as GOES-8, and in 
some channels are slightly superior in quality. 
As of the time of this writing, the plan is to 
move GOES-9 to 135W after a successful 
check-out, where it will provide coverage of 
eastern and central Pacific hurricanes on par 
with the quality of GOES-8. 
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I. Introduction 

In the past, cloud-motion vectors (CMVs) were 
produced within NESDIS in a highly interactive way 
requiring a large investment in manpower. Operators 
displayed loops of geostationary imagery within several 
sectors of the full disk, selected suitable targets on the 
imagery, and followed them manually through 
subsequent images. Heights and vector displacements 
were computed automatically and displayed graphically 
over the imagery for inspection, but height assignment 
methodology was inaccurate when applied to semi- 
transparent cirrus. Operators often made judgments 
on the height of tracers based on their experience, and 
manual cloud heights were often substituted into the 
final product. 

In 1992, the Cooperative Institute for 
Meteorological Satellite Studies (CIMSS) delivered the 
first version of the automated CMV software package to 
NESDIS. It was now possible to produce a full disk 
wind set without manual intervention. Suitable tracers 
were automatically selected within the first image of the 
loop and heights were assigned using several methods. 
Implementation of the C02-slicing algorithm for 
assigning heights to semi-transparent features (Menzel 
et at., 1983) enabled more accurate height assignment 
without operator input. The tracking of features 
through the subsequent imagery was automated using 
the auto-correlation technique (Merrill et al., 1991). It 
was no longer necessary for operators to manually loop 
through imagery and 'point and click* on features. 
Finally, an automated quality-control (QC) algorithm, 
the auto-editor (Hayden , 1993), was developed. 
Initially, the auto-editor was supplemented with manual 
editing. Normal operational procedure involved running 
the automated software, plotting the results, manually 
deleting vectors where appropriate, and then manually 
adding vectors to improve coverage. 

Since 1992, CIMSS and NESDIS have made 
several changes to all aspects of the automated CMV 
software. The most significant of these changes 
involves a new tracer selection procedure, the 
introduction of the H20-intercept height assignment 
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technique, and full automation with a new version of 
the auto-editor. This latest version of the software has 
been running on GOES-8 data in Washington D.C. 
since the fall of 1994. Results have shown that, even 
without the benefit of manual QC, the automated 
GOES-8 cloud-drift winds are superior to any previous 
NESDIS CMV product. 

II. Tracer Selection 

The initial implementation of the tracer 
selection algorithm was very basic. For each target 
domain, the highest pixel brightness values were found, 
and local gradients were computed around those 
locations. Any of those gradients that were large 
enough and not too close to one another were assigned 
as target locations. In practice, the necessary gradient 
threshold was almost always met and the targeting grid 
was very coarse. The result was an extremely 
widespread and uniform field of targets (figure 1 ). 



Figure 1. Typical target distribution resulting from 
GOES-7 imagery and the old tracer selection 
algorithm 


In the new tracer selection algorithm, 
gradients are computed for each pixel within the 
targeting box and the maximum gradients undergo a 
pixel brightness check to ascertain whether they should 


be considered further. The pixel brightness check is 
more a means of filtering out clouds of certain levels 
than quality control. 

All prospective targets now also undergo a 
spatial-coherence analysis (Coakley and Bretherton, 
1982). Data within the target domain is divided into as 
many 3x3 boxes as possible and means and standard 
deviations of the 9 points within these boxes are 
computed. Boxes having standard deviations below a 
given threshold are deemed to be ‘coherent’ and are 
input to a 2- dimensional cluster analysis scheme to 
obtain the mean brightness for any coherent clusters 
(i.e. signals). Presently two checks are employed to 
filter out unwanted targets. First, only two coherent 
signals are allowed. The presence of more than 20 
percent of the coherent 3x3 boxes outside of the two 
largest clusters represents a multi-deck scene. Cloud- 
height diagnosis for multi-deck scenes is extremely 
difficult (Menzel et al. , 1983) so they are eliminated. 
Second, scenes are not allowed to be too coherent. If 
more than 80 percent of the total number of 3x3 boxes 
within the target domain meet the standard deviation 
threshold, the scene is deemed to be too coherent and 
it fails to be a target. Scenes which are that coherent 
typically represent uniform coastal features. 

As a result of the spatial-coherence filtering, 
the total number of targeting attempts within a given 
image can be increased significantly without any 
increase in the overall time necessary to process 
winds. The extra time used to complete a spatial- 
coherence analysis for ail prospective targets is offset 
by the time saved by not having to height assign and 
track undesirable targets. The field of resulting targets 
from the new scheme (figure 2) shows a higher density 
of tracers in desirable locations, and almost none in the 
large clear areas that the previous scheme had 
targeted. 



Figure 2. Typical target distribution resulting from 
GOES-8 imagery and the new tracer selection 
algorithm 


III. H20-lntercept Height Assignment 

As mentioned, the development of the C02- 
slicing algorithm was a major step in the automation of 
wind production. This is the preferred method for 
assigning heights to semi-transparent clouds, but the 
lack of a C02-absorption channel on the current series 
of GOES imagers has necessitated the implementation 
of the H20-intercept algorithm. Comparisons of the 
two methods have demonstrated that the H20-intercept 
algorithm is an adequate replacement (Nieman et. al, 
1993). 

The algorithm is predicated on the fact that the 
radiances for two spectral bands vary linearly with 
cloud amount. Radiances from the infrared window 
and H20-absorption bands are measured and 
compared to calculations of radiances for both of these 
bands for opaque clouds at varying levels in an 
atmosphere specified by a numerical prediction of 
temperature and humidity profiles. Measured and 
calculated radiances will agree for clear-sky and 
opaque cloud conditions. The cloud-top height is 
inferred from the linear extrapolation of measured 
radiances onto the opaque cloud calculations. 

The measured radiances used to infer the 
linear relationship between the two bands are the 
average radiances for the cluster of clearest (warmest) 
fields of view and the cluster of cloudiest (coldest) 
fields of view within the observational area. If the 
calculated water-vapor radiance for clear sky is less 
than the measured water-vapor radiance, the 
calculated water-vapor radiances are adjusted to agree 
with the measured and the difference is attributed to an 
inaccurate transmittance used in the computation of 
dear-column radiance. Calculated warm radiances 
that are greater than measured are not adjusted since 
the low measurement may be the result of cloud 
contamination. 

IV. Automated Quality Control 

An important improvement to the automated 
wind production system came with the new version of 
the automated quality control algorithm (Hayden, 1993) 
shown in figure 3. The basic tool of the objective editor 
remains the 3-dimensional recursive filter objective 
analysis developed at CIMSS. This system has 
undergone an extensive upgrade with attention to 
improving quality control, especially the recursive-filter- 
flag (RFF), which is appended to each datum following 
the analysis. Differences in vector root-mean square 
error (VRMS) with respect to rawinsonde observations 
throughout the range of acceptable RFF flags 
amounted to 0.4 ms' 1 in the old algorithm. The new 
configuration shows an expanded range of 1.1 ms' 1 
from the lowest acceptable RFF values to the highest. 
Thus, the resultant quality flag on the final winds 
should now be more meaningful to the user. 



Figure 3. The objective editing system 


Numerous suggestions from the operational 
community have found their way into the latest version. 
The pressure reassignment is now constrained to 150 
hPa by the algorithm. A small percentage of the 
vectors (< 5 percent) are eliminated by this added 
constraint. Also, the revised editor has tightened the 
tropopause check to a threshold lapse rate of less than 
0.5 K per 25 hPa above 300 hPa in hopes that it will be 
less likely to re-assign heights to some stratospheric 
value. 

The revised editor contains an adjustment for 
the well documented slow bias of CMVs. The quarterly 
statistical summaries of the ECMWF routinely show 
this feature, with respect to both wind measurements 
and coincident speeds of the ECMWF forecast. For 
GOES-7 winds produced with the old automated 
system, the bias could reach -5 ms' 1 at higher wind 
speeds. The bias is also present in comparisons with 
the NMC aviation forecasts, which are used in the 
derivation and quality control of GOES CMVs. To 
mitigate this slow bias, the revised editing procedure 
now increments each vector with 7 percent of the speed 
of the forecast, interpolated to its reassigned level, 
provided that the forecast wind speed is greater than 10 
ms' 1 . Philosophically, this is no different from the bias 
corrections routinely applied to radiance measurements 
to achieve agreement with forward calculations. 

V. Results 

The quality of CMVs is traditionally assessed 
through collocation with rawinsonde (RAOB) 
observations. More elaborate verification procedures 
involve the diagnosis of the impact the inclusion of a 
CMV dataset can have on a numerical model forecast. 
RAOB statistics are more useful for a fixed measure of 


product quality over a period of time. Therefore, that is 
what we present here. 

Figure 4 shows verification statistics for 
GOES-7 cloud-drift winds (G7CD), GOES-8 cloud-drift 
winds (G8CD) and GOES-8 water-vapor motion winds 
(G8WV). VRMS and bias statistics are shown. All 
statistics are monthly means resulting from collocation 
to RAOBs within 2 degrees latitude and longitude and 
25 hPa at 0000 UTC and 1200 UTC daily. All CMVs 
are the product of some version of the automated wind 
software. Operational G7CD winds still undergo a 
manual QC phase in which operators delete some 
vectors; the statistics represent conditions before those 
deletions occurred. G8CD and G8WV winds do not 
undergo go any manual QC. 
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Figure 4. Monthly statistics for VRMS (upper 
curves) and bias (lower curves) for 3 NESDIS 
products 


Statistics for the G7CD winds begin in 
November of 1 992, shortly after the initial 
implementation of the automated software. During that 
first winter, VRMS and bias were high, mainly because 
the auto-editor was mistakenly hot being run at 0000 
UTC daily. This was corrected in the early spring of 
1993 and the statistics improved. The winter of 1993- 
94 was much better in terms of VRMS and bias. The 
old auto-editor was still being run with a few minor 
modifications and it seemed to do a reasonable job. 
There was however, still some degradation in the CMVs 
early in the season. The unedited winds (not shown) 
tend to show an increase in VRMS as the mean wind 
speed increases during the winter. With the advent of 
the first auto-editor this was greatly reduced. The 
entire year of 1994 was the best ever for automated 
cloud-drift wind production as VRMS remained 
consistently at or just above 7ms' 1 . In early 






September 1994 the wind software was shifted from the 
IBM 4381 mainframe to IBM rs6000 workstations. This 
did not adversely affect the quality of the product and 
greatly decreased the amount of time necessary to 
complete a wind set. The new auto-editor also became 
operational at this time (although the bias adjustment 
was not yet being performed) and the quality of the 
G7CD product was as good or better than it had ever 
been for the next four months. In late December of 
1994 the auto-editor began to apply a bias adjustment 
of 5%. This change was immediately apparent as the 
bias dropped almost 1 ms 1 and VRMS values 
improved by approximately 0.15 ms' 1 . Experimentation 
with the G8CD product at CIMSS suggested that even 
better results could be attained with a bias adjustment 
of 7% and that change was implemented for the G7CD 
product in the early summer of 1995. The bias fell well 
below 1 ms' 1 at that time and VRMS values were 
below 7 ms' 1 . 

G8CD wind production began in Washington 
in October of 1994. The GOES-8 products are not yet 
operational so they do not benefit from data acquisition 
priority and backup procedures that are routine for 
GOES-7. Thus the mean statistics are for a widely 
varying number of wind sets depending on availability, 
ranging from 60% to 90%. The initial implementation 
of the G8CD product included the new tracer selection 
procedures, the H20-intercept method in place of the 
C02-sticing algorithm and the new auto-editor using a 
bias adjustment of 7%. The VRMS and bias for the 
G8CD product have been routinely superior to those for 
the G7CD product by about 0.5 ms' 1 . 

G8WV wind production began in Washington 
in December of 1994 using the same software that 
produces the G8CD winds. New tracer selection 
procedures are not enabled for the G8WV product 
since clear air targets in areas of relatively weak water- 
vapor gradients are desired. A simple water-vapor 
brightness temperature comparison scheme is used for 
height assignment (the inclusion of window channel 
data during height assignment offers little in clear air 
since the water vapor channel can penetrate only to 
about 500 hPa). The bias for the G8WV winds is 
comparable to the G8CD winds, while the VRMS 
appears to be consistently greater by 0.5 to 1.0 ms' 1 . 
Even so, VRMSs of 7.5ms r are greatly superior to any 
water-vapor wind product from previous geostationary 
satellites, and it is felt that the inclusion of numerous 
vectors in clear air will complement the coverage 
provided by the G8CD product. 


VI. Summary 

The latest version of the automated cloud 
motion vector software has yielded significant 
improvements in the quality of the GOES cloud-drift 
winds produced operationally by NESDIS. Cloud 
motion vectors resulting from the automated system 
are now equal or superior in quality to those which had 
the benefit of manual quality control a few years ago. 
The single most important factor in this improvement 
has been the upgraded auto-editor, improved tracer 
selection procedures eliminate targets in difficult 
regions and allow a higher target density and therefore 
enhanced coverage in areas of interest. The 
incorporation of the H20-intercept height assignment 
method allows an adequate representation of the 
heights of semi-transparent clouds in the absence of a 
C02-absorption channel. Finally, GOES-8 water-vapor 
motion winds resulting from the automated system are 
superior to any done previously by NESDIS and should 
now be considered as an operational product. 
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The barotropic model is used to explore the advantages of parallel processing in deterministic forecast- 
ing. We apply this model to the track forecasting of hurricane Elena (1985). 

In this particular application, solutions to systems of elliptic equations are the essence of the computa- 
tional mechanics. One set of equations is associated with the decomposition of the wind into irrotational 
and nondivergent components— this determines the initial nondiver gent state. Another set is associated 
with recovery of the streamfunction from the forecasted vorticity. We demonstrate that direct parallel 
methods based on accelerated block cyclic reduction (BCR) significantly reduce the computational time 
required to solve the elliptic equations germane to this decomposition and forecast problem. A 72-h track 
prediction was made using incremental time steps of 16 min on a network of 3000 grid points nominally 
separated by 100 km. The prediction took 30 sec on the 8-processor Alliant FX/8 computer. This was a 
speed-up of 3.7 when compared to the one-processor version. 

The 72-h prediction of Elena’s track was made as the storm moved toward Florida’s west coast Ap- 
proximately 200 km west of Thmpa Bay, Elena executed a dramatic recurvature that ultimately changed 
its course toward the northwest Although the barotropic track forecast was unable to capture the hur- 
ricane’s tight cycloidal looping maneuver, the subsequent northwesterly movement was accurately fore- 
casted as was the location and timing of landfall dear Mobile Bay. 

•) 

* 

i KEY WORDS: Barotropic model, Poisson’s equation, block cyclic reduction, parallel computing. 

C.R. CATEGORIES: G.1.8, J.2, Ell. 


1. INTRODUCTION 

The use of barotropic models to track hurricanes (typhoons) was initiated in Japan 
during the early 1950s (Sasaki and Miyakoda, 1954). Research continued into the 
1960s in both Japan and the U.S.A., where the success of research models eventu- 
ally led to operational implementation at the National Hurricane Center (NHC) in 
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1968 (Sanders et al., 1975). For a history of barotropic modeling as applied to the 
hurricane tracking problem and commentary on its relevance in current research, 
see DeMaria (1985) and Shapiro and Ooyama (1990). 

The sensitivity of barotropic track forecasting to perturbations in the initial con- 
ditions was convincingly demonstrated in the paper by Sanders and Burpee (1968). 
Determination of the initial state is especially difficult due to the scarcity of data 
over the oceans where tropical storms develop. Most observations over the ocean 
are obtained from weather satellites. Typically, these observations are wind esti- 
mates made from a sequence of photographs or images — so-called “cloud-tracked 
winds” or “water-vapor winds” (Velden et al., 1992). Since the barotropic model 
assumes a nondivergent wind, the observations of wind must be decomposed or 
partitioned into nondivergent and irrotational components. 

In an effort to improve the initial conditions for barotropic forecasting, recent 
efforts have been made to include a time history of wind analysis into the estimate 
of the initial state. To accomplish this task in a dynamically consistent fashion, ad- 
joint methods have been used (Lewis et al., 1987; DeMaria and Jones, 1991). These 
assimilation/initialization methods require a model of comparable complexity to the 
forecast model. The so-called adjoint model is integrated backward in time and is 
used iteratively with the forward model (forecast model) to find an optimal initial 
condition that minimizes the discrepancy between model forecast and observations. 
Although aesthetically pleasing and theoretically attractive, the computational de- 
mands have made this strategy prohibitive in the operational environment With the 
research models tested to date, the computational cost of an adjoint assimilation 
with barotropic models is « 10-20 times as costly as the forecast alone (Lewis and 
Derber, 1985; Talagrand and Courtier, 1987). Realistic implementation of adjoint 
assimilation in operations will require the maturing technology of multicomputer 
and multiprocessors (Lakshmivarahan and Dhall, 1990). We refer to this technology 
as parallel processing or parallel computing. 

Our ultimate aim is to apply parallel computing technology to both the forecast- 
ing and assimilation problem (adjoint assimilation). We concentrate on the forward 
problem in this initial effort, but the transfer of the computational algorithms to 
the backward or adjoint problem should be straightforward. We apply our parallel 
processing mechanics to the track forecast of hurricane Elena (1985). This case is 
especially challenging because of the looping maneuver that the storm executed in 
the northeastern Gulf of Mexico. This case study also benefits from an unusually 
large set of synoptic-scale wind analyses produced at the Space Science and Engi- 
neering Center (University of Wisconsin-Madison). 

In Section 2, following Lynch (1988, 1989), we describe a computational scheme 
for decomposing a windfield into non-divergent and irrotational components. This 
scheme leads to the solution of two related Poisson’s equations. The baratropic 
model itself is described in Section 3. This has two components— one is the time- 
dependent evolution of vorticity and the second is a Poisson’s equation relating the 
stream function and vorticity. Since the solution of Poisson’s equation is central to 
the computations, we use a direct parallel method, called the block cyclic reduction, 
for solving these equations. A summary of this method is given in Appendix A Sec- 
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tion 4 discusses various issues related to choice of grid, boundary conditions, nature 
and type of discrete approximations, and conditions for stability. Speed-up analy- 
sis resulting from solving this model on the state of the art multivector processor 
Alliant FX/8 is given in Section 5. Concluding statements are given in Section 6. 


2. DECOMPOSITION 

To facilitate our work, we have mapped the earth onto the Mercator projection. 
In the kinematic and dynamics equations that follow, the east-west and north- 
south coordinates — x and y, respectively — are distances on the Mercator plane and 
the image scale factor, m, is the ratio of distance on this projection to the corre- 
sponding distance on earth. This image factor is given by (cas<fo/ cos<p) where <po is 
the latitude of intersection between the Mercator cylinder and the spherical earth 
( 0o = 30° N in our case). Saucier (1955) discusses the geometry of the Mercator 
projection and includes the derivation of the image scale factor. 

Let V denote the horizontal wind vector, x the velocity potential, $ the stream 
function, ( the vorticity and 6 the divergence. Let V^, and \ x be the nondivergent 
and the irrotational components of this wind field, respectively. Helmholtz’s theo- 


rem permits decomposition of V as follows: 


V = V^, + V*, 

(1) 

where 


V x = k x Vtf, 

(2) 

ll 

<1 

X 

(3) 

and V = ( m(dldx),m(d{dy )) is the 2 — d gradient operator. To simplify notation, 
the image-scale factor will be assumed to be associated with all differential opera- 
tors, but it will not be explicitly written. Thus, d/dx is in actuality m(d/dx), (d^/dx 2 ) 
means m(d/dx)(m(d/dx)), etc. Let i, j and k be unit vectors in the 2 orthog- 
onal horizontal directions (x,y) and the vertical direction (z), respectively. Since 
V • 0 and V x V x = 0, it follows that 

( = lcVxV = k-Vx(kxV$): 

= V 2 * (4) 

and 

S = V-V = V 2 x 

(5) 

where 

v * = (iL + iL) 

\a* 2 fly 2 / 


Since V = in + jv, then 



( =k.v* v =(!-D 


( 6 ) 
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Table 1 Lynch’s Eight Different Variations of Boundary Conditions for Decomposing the Wind Field 
(Dirichlet and Neuman conditions are denoted by D and N, respectively) 


Version 

Number 

First Boundary 
Condition 

Second Boundary 
Condition 

Type of First 
and Second B.G 

1 

X 

II 

o 


DD 

2 

dn 7 " 


ND 

3 

* = 0 

m ^ - V 

m Tn~ Vn 

DN 

4 

8|* 

II 

* 

nt — v n + m — 

dn ds 

NN 

5 

x = 0 

II 

DN 

6 

81# 

II 

8|S 

II 

1 

$\8 

NN 

7 

¥ =0 

j;) ds 

DD 

8 

?!* 

II 

3 

x = ffa-ys)ds 

ND 

and 


■ ( du dv\ 

6 ~\d^ + dy)' 

(7) 


To extract the nondivergent wind component from the observed wind, we first 
compute the vorticity and divergence which serve as the forcing functions in the so- 
lution to Poisson’s equations, Eqs. (4) and (5). These equations are coupled through 
the boundary condition (b.c.’s) which can assume a variety of forms. As discussed 
by Lynch and collaborators (Bijlsma et al., 1986; Lynch, 1988, 1989), the decompo- 
sition is not unique in a limited domain and a wide variety of results are possible 
that depend on the choice of b.c.’s. Table 1 displays eight different sets of b.c.’s that 
have been tested by Lynch (1988). The first four versions use the normal compo- 
nent (V n h) of the observed wind on the boundary whereas the last four use the 
tangential component (Fji). The parameters 

7 "( = § s r/ v "‘ <s ) “ d »(-£-r/ v '") 

are constant gradient conditions on the boundary that satisfy divergence theorem 
identities — so-called compatibility conditions. Section 2 of Lynch (1988) can be read 
for a more detailed discussion of these boundary condition choices. 

Our choices are narrowed by choosing b.c.’s that minimiz e the kinetic energy in 
the divergent component ( K x ). This is dictated by the desire to minimize the kinetic 
energy in the “largely spurious divergent part of the wind” (Sanders and Burpee, 
1968). As shown by Davies-Jones (1988) and Lynch (1988), x = 0 on the boundary 
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is sufficient to minimize K x and furthermore causes the nonphysical cross-term K$ x 
to vanish (again see Lynch (1988) for details). Version 1 (due to Sangster (I960)) 
and version 5 are, therefore, candidates. Dirichlet b.c.’s are easier to implement 
than Neumann conditions since the solution is unique, i.e., it does not involve an 
arbitrary additive constant, nor is there a compatibility condition. Thus, version 1 is 
preferred. Care must be exercised, however, to insure that 


is single valued; subscript B is used to indicate the value on the boundary. This 
requirement can be satisfied by first solving (5) for x w ith Xb - 0, computing 
dx/dn | B and adjusting V„ (to V„) by a constant amount e around the boundary 
such that 


L(rn- v '-‘) ds -°- 


The final step is to solve (4) with 



where 'I'b(So) = 0 


by choice. Thus, version 1 has the same desirable properties as version 5 without 
the complexities of a Neumann b.c. The only difference is that version 1 uses V„ 
on the boundary while version 5 uses Vs. There doesn’t appear to be any physical 
reason for preferring V$ over V„, or vice versa. 


3. THE BAROTROPIC MODEL 

The approach to barotropic forecasting follows Leith (1978) and DeMaria (1987). 
The prognostic equation is given by 

§-j«+/,n (8) 

where J is the Jacobian and / is the Coriolis parameter. 

In this prognostic equation, ( represents a modified vorticity given by 

£ = (V 2 -<t 2 )*. (9) 

The derivation of (8) is found in Leith (1978) and follows from operations on the 
shallow water equations. As discussed in DeMaria (1987), this form of the equation 
amounts to a divergence-corrected barotropic model which prevents the retrogres- 
sion of ultralong Rossby waves. The parameter a is empirically determined and 
is taken to be l/800(km -1 ) based on a recommendation from DeMaria (personal 
communication). We still use the nondivergent wind to initialize the model. 

Given the wind components (u, v) at each of the grid points, a scheme for com- 
puting the solution of (8) is described as follows: 
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(a) Compute ( and 6 from the initial wind field, then solve (4) and (5) for and 
X following the decomposition process described in Section 2. 

(b) Use (8) to step forward and find ( at the next incremental time. The vorticity 
at the boundary is calculated differently, depending on whether there is inflow 
or outflow at the boundary. If outflow, the vorticity is calculated by a one- 
sided finite difference version of (8). If inflow, the vorticity is unchanged. 

(c) Using the forecasted vorticity, (9) is solved for 

(d) Repeat steps (b) and (c) until the entire forecast period is covered. 

Both the decomposition of the wind field and the barotropic forecast hinge on 
the solution to Poisson’s equation, lb find the solution to this equation, we apply 
a direct method called block cyclic reduction (BCR). Since the use of this method 
in a parallel processing mode is critical to our work, salient features of this method 
are given in the appendix. For more details, refer to Buzbee, Golub and Nielson 
(1970), Hockney (1965), Lakshmivarahan and Dhall (1991), Swartztrauber (1984), 
and Sweet (1977, 1988). 


4. APPLICATION TO HURRICANE ELENA 

Hurricane Elena (1985) had its origin in a rapidly moving tropical wave that was 
first analyzed over the Saharan Desert (Velden, 1987). It moved over the Atlantic 
Ocean at a rate of « 15 ms -1 and appeared as a tropical depression over Cuba on 
August 28th. The depression intensified and was classified as a tropical storm by 
00 UTC on the 29th. It was upgraded to hurricane status 12 hours later when it 
was centered at ~ 25N, 85W. Subsequent positions at 12 h intervals are indicated by 
dots along the solid line in Figure 1. The initial time for our forecast was 12 UTC 
31 August, when the storm occupied the position indicated by (§ ) in Figure 1. 

a. Numerics: All experiments were conducted on the ALLIANT FX/8 using 
double precision arithmetic. Five separate modes of computation were used starting 
with the scalar mode — use of a single processor without the benefit of vectorization. 
This was followed by computation with 1, 2, 4, and 8 processors with vectorization. 
Since the BCR-partial fraction method was more efficient than the polynomial fac- 
torization form of BCR (Jwo et al., 1992), all experiments used the partial fraction 
approach (see Appendix). 

The equations are solved on the staggered mesh of grid points displayed in Fig- 
ure 2. This staggering is in accord with Arakawa’s C-grid convention (Arakawa, 
1966). The finite-difference approximation to the Jacobian conserves the vorticity, 
the kinetic energy, and the squared vorticity. 

Following DeMaria (1985, 1987), we adopted the Adams-Bashford scheme for 
integrating the model (Anderson et al., 1984). As suggested by DeMaria (personal 
communication), a reasonable limitation on the time-step is given by 

A < A 

‘ ~ 2m\/2(M 2 + v 2 )^ 


(25 min), 
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Figure X The observed and forecasted tracks of Elena (1985). The location of Elena at the initial time 
(t = 0, 12 UTC Aug. 31) for barotropic forecasts is indicated by ( ). The numbers 1, 2, and 3 indicate the 
location of Elena at 1, 2, and 3 days after t = 0. Our forecast is termed barotropic and the operational 
model is SANBAR. 



li 

Vi 

* 


* 


* 


b X* 1 

Figure 2 The lower right hand corner of the staggered grid used in our study. The boundary values are 
denoted by the use of subscript 1)’. 


where max signifies the largest wind speed, A and A t are the space and time incre- 
ments, and m is the image scale factor at the location of the maximum speed. As 
indicated in parentheses, the limiting A, is 25 min for our grid and wind field. A 
time step of 16 minutes was used in all of our computations. 
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Prior to the forecast, an artificial axisymmetric vortex — representing the hurri- 
cane circulation — is superimposed on the 'P-field found from the decomposition 
process. This helps us locate the hurricane at subsequent times. The vortex remains 
intact through the 72 h forecast period, but there is a modest amount of numerical 
smoothing that weakens the vortex. The mathematical form of the vortex is: 

V„ W = V„(l-)exp{i[l-(^)]} and faW-; 

( 10 ) 

where V w is the tangential wind and (# the associated vorticity. The parameters in 
the wind equation are: V m (the maximum tangential wind), r (the distance from the 
storm center), r m (the radius of maximum wind), and b (an empirical parameter that 
determines the rate at which V# decreases at large radii). We have used DeMaria’s 
values for these parameters to typify observed Atlantic tropical storms (DeMaria, 
1987). These values are V m - 25 ms -1 and r m = 150 km, and the parameter b is 
chosen such that V# = 5 ms -1 at 600 km (this gives a value of b = 0.998). 

b. TVuck forecast: Barotropic forecasting for a particular level in the tropical 
atmosphere has been unsuccessful for a variety of reasons that have been enumer- 
ated in Sanders and Burpee (1968). They determined “ . . . that forecasting should 
be carried out for the tropospheric mean flow ... rather than by selection of a 
single, presumably representative level.” This work was the foundation of the op- 
erational model developed at NHC (Sanders et al., 1975). The operational model 
has become known as the Sanders Barotropic model, or SANBAR in acronymic 
form. Following the work of Sanders and Burpee (1968) and Sanders et al. (1975), 
we calculate a mean wind throughout the troposphere using satellite-derived winds 
and radiosondes following the procedures established at the Space Science and En- 
gineering Center in Madison, Wisconsin (UW-SSEC) (Velden et al„ 1984). In sub- 
sequent discussion, we refer to these winds as deep-layer-mean (DLM) winds. 

Figure 3 shows the DLM wind field at t = 0 (12 UTC 31 Aug.) plotted over 
our analysis and forecast domain. A well-defined trough that extends from New 
England down through the northeast comer of the Gulf of Mexico is apparent in 
the DLM wind vectors. The initial streamfunction field associated with the DLM 
wind is shown in Figure 4. The artificial hurricane vortex described in Section 4a 
has been superimposed on this large-scale $ field. The forecasted stream function 
fields at t = 72 h is shown in Figure 5. The path of the hurricane over the three-day 
period has been superimposed on the 72 h forecast of streamfunction. 

We estimated Elena's position in our forecasts by finding the center of the imbed- 
ded vortex at 12 h intervals. This is done by finding the minimum value or the 
maximum £ value. We found insignificant difference in positioning Elena by ei- 
ther method and the forecasted positions at 12 h intervals are displayed in Fig- 
ure 1. The observed track as well as the forecasted track from SANBAR are also 
shown in this figure. The observed tight cycloidal motion (re-curvature) between 
t — 0 and t = 24 h was not captured by our model forecast. This small scale feature 
in the trajectory is near the limit of the model resolution, i.e., the area outlined 
by this “teardrop” shaped part of the trajectory is « (0.5° lat) 2 = (50 km) 2 . A grid 
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Figure 3 Distribution of deep-layer-mean winds (DLM) at t = 0 (12 UTC, 31 Aug 85). The wind speed 
is scaled by the largest speed which is shown in the lower right hand comer (zz 30 ms -1 ). Elena’s posi- 
tion at this time is indicated by ( ). 


c I»Ji 



Figure 4 Streamfiinction associated with the DLM winds and imbedded vortex at t = 0. The contours 
of the streamfunction are drawn at intervals of 10 6 m 2 s” 1 . 
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Figure 5 Same as Figure 4 except map is valid at t = 72 h. The observed track of Elena between 0-72 h 
has been superimposed. 


cell in the model is ss (100 km) 2 . Although the forecasted location of landfall was 
% 150 km east of the observed location (near Mobile Bay), the time for landfall was 
accurate. The positioning error was less than the average error for 48 h predictions 
(NHC, personal communication). 

The parallel-barotropic forecast was an improvement over SANBAR because it 
did not use the “initial storm motion” (ISM) vector as part of the forecast. It has 
been found that operational forecasts with SANBAR generally benefit from the 
inclusion of an ISM vector based on the storm’s previous path. In cases exhibit- 
ing re-curvature, however, incorporation of recent history (persistence) is generally 
detrimental to the track forecast Since re-curvature is difficult to anticipate and 
since it is a relatively rare event, the operational forecasts always include the ISM 
vector in the initial condition. Examination of the forecast from SANBAR indicates 
that the imposition of the ISM constraint at t = 0 forced Elena to erroneously move 
over the northern portion of the Florida peninsula during the early part of the 72 h 
forecast. 

5. SPEED-UP ANALYSIS 

Table 2 summarizes the results for our barotropic forecast using the Alliant multi- 
vector processor. This table illustrates two types of speed-up computations. The 
speed-up ratio 1 is the ratio with respect to the serial time. Using one vector pro- 
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Table 2 Speed-Up Ratios for 72 h Forecasting Process 


Number of Processors 

Elapsed Time (Sec) 

Speed-Up Ratio 1 

Speed-Up Ratio 2 

Serial 

216.18 

1.0 



1 Vector 

109.47 

1.97 

1.0 

2 Vector 

61.03 

3.54 

1.79 

4 Vector 

38.82 

557 

2.82 

8 Vector 

29.53 

732 

3.71 


cessor we obtain a speed-up of nearly 2, vising 8 vector processors a speed-up of 
7.32. In the latter case, there are two factors contributing to the speed-up, namely, 
parallelism (8 processors) and vectorization in each of the eight processors. In con- 
trast, the speed-up ratio 2 measures the speed-up resulting from multivectors to 
single vector processing. Thus, eight vector processors result in a speed-up ratio of 
3.71. Since vectorization is a common factor in the second ratio, the effect of paral- 
lelism on speed-up is isolated. 

Referring to the results in Table 2, one might wonder why a speed-up ratio of 
ss 3.7 is achieved when 8 processors are used. In practice, a number of factors af- 
fect the speed-up achievable. First, the algorithm itself is a factor. It could happen 
that there are segments of the computation that are intrinsically serial and cannot 
be parallelized at all. This will lead to loss of speed-up as dictated by the well- 
known Amdahl’s law (Lakshmivarahan and DhaU, 1990). Referring to the model 
equations (8) and (9), it can be verified that (refer to Appendix A), while BCR 
method for solving (9) exhibits “good” parallelism, the degree of parallelism does 
not remain constant, which is an intrinsic property of this class of algorithms. In- 
deed, the degree of parallelism decreases during the course of the reduction phase, 
and then increases in the back-substitution phase. In other words, the processor uti- 
lization decreases and increases during the reduction and back-substitution phases 
respectively. Secondly, stability of numerical integration dictates the use of Adams- 
Bashford scheme in combination with Arakawa’s approximation for the Jacobian. 
However, the computational process resulting from the use of these schemes do not 
readily vectorize and/or parallelize. These factors limit the speed-up that is achiev- 
able. 

Second, the machine architecture is a factor. In the Alliant FX/8, for example, 
the processors are connected to the main memory through a fast cache memory. 
The size of this cache limits the data transfer rate between the memory and the 
“thirsty” processors. Thus, when all eight processors are concurrently working in 
the vector mode, the data needed by these processors may not be readily available 
in cache memory which may affect the data transfer rate between memory and the 
processors. Thus, in such a situation, increase in number of processors results in 
diminishing returns in terms of the speed-up. Increasing the cache size will only 
postpone this saturation effect, but on architectures with no cache memory, such as 
the CRAY research machines, linear increase in speed-up can be achieved. 

Thirdly, when one increases the number of processors, unless the problem size is 
also increased, the given quantum of work is merely smeared across a large num- 
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ber of processors. Oftentimes it may take more time to distribute the data to the 
processors than to compute the results. This contributes to idling of processors and 
in turn reduces the processor efficiency, which is speed-up per processor. Since the 
data set available on Elena was fixed a priori, increasing the number of processors 
would only result in diminishing returns as evidenced in Thble 2. Thus, success in 
parallel computing critically depends on the proper match between the algorithm, 
the architecture, and the problem size. 

6. CONCLUSION 

We have taken the simplest deterministic weather prediction model that is still used 
operationally — the barotropic model — and have brought the maturing technology 
of parallel processing to bear on its solution. We have applied the computational 
mechanics to the problem of tracking hurricane Elena (1985). 

The results are encouraging from both a pragmatic computational viewpoint as 
well as a practical weather forecasting viewpoint. A significant saving of compu- 
tational resources is possible with parallel processing. Essential to the time-saving 
is the use of parallel direct methods for solving the elliptic equations associated 
with the barotropic forecast We have found that the block cyclic reduction (BCR) 
method using partial fraction expansion is an efficient strategy. As a measure of the 
efficiency, we calculated the speed-up ratios with respect to both the serial mode 
and the one vector processor mode and found that a speed-up ratio of 3-4 is possi- 
ble when an 8 vector processor is compared to a single vector machine. 

A complete set of DLM winds for Elena (1985) were derived at the UW-SSEC 
and used as input to the barotropic model. A decomposition of these winds was 
accomplished in such a way that the kinetic energy in the nondivergent component 
was maximized. The streamfunction extracted from the decomposition process was 
used to initialize the barotropic model and application was made to the period of 
time when Elena executed a dramatic re-curvature in the eastern Gulf of Mexico. 
Although the model did not capture the small scale cycloidal motion associated 
with the abrupt change in Elena’s path, the larger-scale redirection from easterly 
to northwesterly was well handled. Especially encouraging was an accurate 48 h 
prediction of landfall. 

Based on this success, we are prepared to address the problem of coding the 
adjoint model that can be used in model sensitivity studies and data assimilation. 
Since the adjoint equations also hinge on the solution to Poisson’s equation, the 
implementation should proceed expeditiously. 

Appendix. Solution to Poisson’s Equation Using a Parallel Processing Version of 
Block Cyclic Reduction (BCR) 

BCR methods were first discussed by Hockney (1965) and were further developed 
by Buzbee et al. (1970). However, applicability to large-scale problems was ham- 
pered by the presence of serial bottleneck which inhibited the exploitation of full 
parallelism (Lakshmivarahan and Dhall, 1990). Only recently have methods been 
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developed to overcome this bottleneck (Sweet, 1988; Gallopoulos and Saad, 1989). 
The classical BCR, when combined with acceleration schemes, gives rise to an algo- 
rithm which holds promise for both parallelism and vectorization (Jwo et ah, 1992). 
We have now developed codes for the BCR method that include two basic options: 
polynomial factorization and the partial fraction expansion. Both versions can ac- 
commodate three types of boundary conditions: Dirichlet, Neumann, or periodic. 
They can also handle an arbitrary M xN grid as well as special grids with M 
(or N) = 2 k — 1, k an integer (Sweet, 1977; Jwo et al., 1992). Poisson’s equation 
is discretized on a uniform Mercator grid (100 km spacing at latitude 30° N). In the 
experiments, variables on the staggered grid have the following dimensions ( M,N ): 
u(65,43), v(64,44), $(65,44) and *(66, 45) where M is the index in the east-west 
direction and N is the corresponding index in the north-south direction. 

Poisson’s equation on an M x N grid is represented by 


Ax = b, 


(A.1) 


where A is a sparse symmetric matrix. When Dirichlet boundary conditions are used, 
A assumes the form: 
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(A.2) 


where there are ( N - 2) block matrices B along the principal diagonal of A, and B 
is a tridiagonal matrix of order (M - 2) x (M - 2). The form of B is: 
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(A.3) 


and C is the identity matrix of order (Af - 2) x (M - 2). 

We shall illustrate block cyclic reduction in the case where A is a 7 x 7 block 
tridiagonal matrix. The system Ax = b can be rewritten as 


Bxi + X2 = bi 


x,_i + Bx, + x,+i = b, 2 < * < 6 


(A.4) 


X6+X7 =by. 
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Consider the sets of three consecutive equations centered around x,, for * = 2, 4 
and 6 

Xj _2 + Bx,_ix, = bj_i ] 


x,_iBxf + %i +1 = b i 
Xi + Bxi+i + Xj+2 = b,+i 


i = 2,4,6, 


(A.5) 


where xo = x% - 0 . 

Multiplying the second equation by — B and adding the three equations of the 
set, we eliminate the odd indexed terms x,_i and x,+i. These operations can be 
accomplished in parallel for i = 2, 4 and 6 , and we obtain 


Xj_2 + (21 - B 2 )x, + Xj+2 = b,+i + b,_i - Bbj, i = 2, 4, 6. 


(A. 6 ) 


This constitutes the first step of the reduction process and results in the split of 
the 7 x 7 system into two subsystems — one for odd indexed Xj’s called the eliminated 
system and the other for even indexed x, ’s called the reduced system 

B(1)X 2 + X 4 = b 2 

X 2 + B(1)X4 + X 6 = b 4 (A.7) 

X4 + B(l)x$ = b^ 

is the reduced system where B(l) = 21 - B(0) 2 , (B(0) = B), which is a matrix poly- 
nomial in B of degree 2 and b,(l) = b,+i + b,_i - B(0)b,, for i = 2, 4 and 6. The 
eliminated system is given by 


B(0) 0 

0 B(0) 

0 0 


bi-x 2 

b3 - X2 - X4 
bs — X4 - X6 
-by-xe 


(A. 8 ) 


Thus solving (A.7) for b 2 , b 4 , and be, we can readily solve for xi, x 3 , x 5 and X 7 from 
the eliminated system. Notice (A.7) is again a block tridiagonal system and we can 
once again multiply the second equation in (A.7) by -B(l) and add the equations 
together to eliminate X 2 and xe. The new reduced system is 


B(2)X4 = b 4 (2), 


(A.9) 


where 


B(2) = 2I-[B(1)] 2 , 

which is a matrix polynomial in B of degree 4 and 

b 4 (2) = b 2 (l) + b 6 (l)-B(l)b 4 (l). 
The new eliminated system is 

B(l) 0 1 rx 2 l = rb2(l)-*4‘ 

. 0 B(1)J [x^ J [b 6 (l)-x 4 . 


(A. 10) 
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1 2 3 4 5 6 7 



Figure 6 Pictorial display of the block cyclic reduction (BCR) process for a linear system of 7 equations 
where each node in the diagram represents an equation in one of the variables: X\,X 2 ,...,X 2 . 

This algorithm is diagrammatically illustrated in Figure 6. It consists of two 
phases. In the first step of the reduction phase, eliminate xi, X3, x$, and X7 in parallel 
to get a reduced system in X2, X4, and X6. In the second step eliminate X2 and x$ in 
parallel to obtain the reduced system in X4. By solving this system we can recover 
other unknowns in the back-substitution phase. Knowing X4, X2 and x$ are recovered 
first in parallel using (A.10) and then from (A.8), xi, X3, X5, and X7 are recovered in 
parallel. 

Computationally, the back substitution involves solving 


B(2)x4=b 4 (2) 

(A.11) 

B(l)x, =b,(l)-X4, for 1 = 2,6, 

(A. 12) 

B(0)x, = b, — Xj+i - x,_i 

(A- 13) 


for i = 1 , 3, 5, and 7, where xo = xg = 0. 

Consider the system in (A.11). From B(2) = 21 - [B(l)] z = 21 - (21 — [B(0)] 2 ) 2 , it 
follows that computation of B(2) involves expensive matrix multiplications. Besides, 
while B(0) is tridiagonal, [B(0)] 2 is not tridiagonal. Thus B(l) and B(2) do not in- 
herit the sparse structure of B(0). Thus solving (A.11) by directly computing B(2) 
is not computationally worthwhile. An alternate approach is to exploit the structure 
of the matrix polynomial and to cleverly factor it so that each factor inherits the 
tridiagonal structure of B which can then be exploited in computation. To this end, 
define a class of polynomials recursively as follows: 

p 2 r+i(a,t) = 2 r 2 '* 1 - [p2r(a,t )] 2 (A. 14) 

where 

p 2 (a,t) = 2t z - a 2 . 

p 2 r(a,l) is known as the Chebyshev polynomial of the first kind and degree 2 r . The 
importance of this class of polynomials stems from the fact that 


®( r ) = P2'(o,t ) \a=B,t=l • 
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UMS, B(l) = 21 - B 2 . By setting a = -2tcos8, it can be shown that 


Clearly, 


P 2 '(a,t ) = 2t r cos 2 r 6. 


P2'(a,t) = 0, 


£“- cos (^rr) ,r ' / = 1,2,-, 2'- 


Using this we obtain a factorization 


= - IK* + 2/ + cos0 y (r)), 


where 


W- (^r)«- 


The corresponding factorization for B(r) is given by 

v 

b(o=-II h ,M' 

;=i 


where 


H ; (r) = (B + 21 cos 0;(r)). 


Notice Hy(r) is tridiagonal since B is tridiagonal and 2Icos0/(r) is diagonal. 

Given this factorization, there are two possibilities for solving systems of the type 

B(r)y = d for r = 0,1,2, We illustrate this using (A.U). 

Recall 

B(2) = -Hi(2)H 2 (2)H 3 (2)H4(2). 

Let Zo = b 4 ( 2 ). Then (A.ll) becomes 

Hi(2)H 2 (2)H 3 (2)H4(2)x 4 = Zo. 

Solve 

H,(2)Z, = Z,_i 

successively for i = 1, 2, 3 and 4. Clearly X 4 = Z r . The advantage of this approach 
called the polynomial factorization is that each H, (2) is sparse. Yet the collection of 
systems is solved sequentially which inhibits parallelism. 

In search of a method that admits more parallelism, rewrite (A. 11 ) as 

*4 = [B(2)]- 1 b 4 (2) 

Expressing [p 2 '(a,0 ] -1 in partial fraction expansion we obtain 

X 4 = [a 1 Hf 1 (2) + a 2 H 2 - 1 (2) + a 3 H 3 " 1 (2) + a 4 H 4 - 1 (2)]b 4 (2). 
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In this expression each of the terms on the right hand side can be independently 
computed in parallel. Swarztrauber (1984) and Lakshmivarahan and Dhall (1990) 
can be consulted for further elaboration on the steps involved in parallelization of 
this algorithm. 
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ABSTRACT 

The first six months of operation with the new automatic cloud motion wind 
(CMW) processing at NESDIS is evaluated. Improved vector accuracy is shown, 
but a slow bias to upper level winds persists. Research development of 
infrared window and water vapor winds using the METEOSAT-3 are discussed. 
With the METEOSAT data, C02 height assignment cannot be used, and the 
infrared window CMW are slightly inferior to those from GOBS- 7. METEOSAT- 3 
water vapor winds are quite encouraging and appear to offer the long- sought 
useful information at middle levels, although the quality is marginal for 
improving on the NMC forecasts. Operational gradient winds generated from 
VAS temperature profiles are also explored, and shown to offer positive 
information on the geopotential gradients. 

INTRODUCTION 

Cloud Motion Winds (CMW) continue to be an important topic of research at 
centers which process geostationary satellite data. In the past several 
years there has been a renewed interest in international cooperation 
resulting in a Scientific Meeting held by COSPAR at the Hague 28-29 June, 
1990; a Workshop on wind extraction from operational meteorological satellite 
data held in Washington D.C., 17-19 September, 1991; a follow-up Workshop on 
the improvement of operational NESDIS CMW 20-21 February, 1992; and now the 
COSPAR meeting in Washington D.C., 4 September 1992. Impetus for the 
increased interest has been generated primarily by the European Center for 
Medium- range Weather Forecasts (ECMWF) and the National Meteorological Center 
(NMC) who impose increasingly stringent requirements for wind observations as 
the sophistication and accuracy of their forecast models improve. It is well 
known that this evolution has outpaced the improvements in accuracy attained 
by satellite data /l/ to the point where the data are, in some cases, no 
longer used. Naturally this is a challenge to those generating the data. 
Beyond this, however, it is also well know that the data from the various 
producing centers have not been of commensurate quality, suggesting that a 
pooling of knowledge might benefit all. The last factor is the main 
motivation for the continuation of the international scientific study groups 
and workshops. 

In this paper we shall address the recent research efforts contributing to 
the operational CMW production at the NESDIS. This research mainly 
accomplished at the Cooperative Institute for Meteorological Satellite 
Studies (CIMSS) and at the Cooperative Institute for Research in the 
Atmosphere (CIRA) where NESDIS maintains groups associated with the 
Universities of Wisconsin and Colorado State respectively. The two 
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institutes have different approaches to the task, CIMSS being primarily 
involved with synoptic flow regimes whereas CIRA concentrates on mesoscale 
applications. Specific topics treated here include: an assessment of the 
recent (February 92) implementation of an automatic CMW production at NESDIS; 
current research in the generation of CMW from METEOSAT (in cooperation with 
ESOC) ; continuing investigation of the utility of winds from water vapor 
imagery; and, finally, recent application of gradient winds derived from VAS 
retrievals of temperature. The last is, to be sure, not a wind derived from 
amimated imagery, but it shares the common bond of being a satellite product 
which is available with high spatial density. Also, the gradient winds have 
been of some interest to the NMC as a compliment (in cloud- free areas) to the 
CMW. 


NESDIS AUTOMATED CMW PRODUCTION FROM GOES 

The automated CMW procedures implemented at NESDIS in February, 1992 
include automatic target selection, initial tracer height assignment by a C02 
slicing algorithm (or some variamt of am infrared window measurement matched 
to a T (p) profile from the NMC forecast if this results in a lower pressure), 
amd assimilation/quality control which may reassign the height by 
assimilation with other vectors amd the NMC 12 -hour forecast. These 
procedures are described in Merrill et al. /2/, Hayden and Velden /3/ amd 
also in the Proceedings of the Workshop on Wind Extraction from Operational 
Meteorological Satellite Data /4/ (hereafter referred to as Workshop Report) . 
The implementation, replacing manual tracking on a sequence of images, 
resulted in a large increase in the number of CMW and an improvement in 
timeliness (imagery from 30 to 90 minutes before analysis time instead of 90 
to 150 minutes before) . Accuracy statistics (as compared with co- located 
rawinsondes) have been collected over the first 6 months of automated 
operation and are presented in Table 1. The Table includes statistics for 
the fully automated sample and also for the "final" sample after manual 
editing. Table 2 shows statistics collected over a comparable period a year 
earlier for the manual CMW. Comparison of Tables 1 and 2 reveals that the 
vector error for high level winds has been reduced (in the final product) 
with the new system from 8.6 to 7.5 ms-1. The improvement is achieved even 
though the mean speed of the later sample is greater, 23 vs. 17 ms-1. The 
manual editing step is effective in reducing the error from 7.8 ms-1. Of 
more importance than the actual magnitude of the error is the capacity of the 
CMW to improve over the NMC forecast. This feature is seen for both old and 
new systems with respect to the high level winds. The mid- level are more 
problematical. In the new system the vector rms is not improved over the 
guess, though it is comparable. Also at mid- levels, the manual editing seems 
ineffective since it does not improve the accuracy relative to the forecast 
although it does lower the magnitude of the vector error. We note that the 
sample of mid- level vectors is considerably larger in the automated system, 
suggesting that the manual target selection is more discriminating in this 
difficult layer where large shears predominate. 

The tables also show a difference in forecast quality between the two 
samples. This is probably the result of using different NMC forecast models 
in the respective processing. For the manual derivation the global forecast 
was used, but with the automated version the aviation forecast has been 
substituted. Particularly with respect to mean speed, the aviation forecast 
appears to be a considerable improvement, and it is encouraging that the CMW 
can improve on this as well. The improved forecast is possibly responsible 
for the failure of automated mid- level vectors to improve on the forecast. 


Deriving Winds from Satellite Data at NESDIS /CIMSS 


(3)101 


TABLES 1-3 

A long standing problem with CMW has been the 'slow bias" which has 
been shown to increase approximately linearly with wind speed as measured by 
comparison to co- located aircraft reports. As shown in the Workshop Report 
/ 4/ , the bias error for the NESDIS manual vector varies from about 5 ms-l for 
the slow category (0-10 ms-l) to about minus 12 ms-1 for the fast category 
(>50 ms-l) . For the automated processing, Fig. 1 presents a scatter plot of 
the CMW speed error (points) and the Forecast error (crosses) vs. rawinsonde 
speed for the manually edited sample of Table 1. Clearly the previously 
noted bias is retained in the automated processing, with a somewhat reduced 
magnitude. This slow bias has also been confirmed using aircraft reports at 
CIMSS, the BCMWF, and the NMC. 

Pig 1 

Fig. i also shows that the NMC forecast, as it is processed in the CMW 
procedures, exhibits a similar, but more serious bias. Our investigation has 
shown that a major part of the apparent forecast bias is caused by a 
reformatting in which forecast grids are interpolated to a data file 
organized as vertical profiles of vectors. (Table 3). However, at the 
higher wind speeds the original fields are also markedly slow, as shown in 
Table 3 . Because the CMW processed forecasts exhibit this bias , the 
assimilation quality control step, which blends the forecast and the CMW, 
does flag the slow CMW. One might even ask if the forecast bias does not 
exacerbate the problem by having heights reassigned to lower wind speed 
regimes. We have examined this latter possibility and determined that it is 
not occurring. Steps are being taken to correct the forecast bias and 
perhaps improve the quality control of the CMW. 

Conclusions which we draw regarding the current NESDIS operational 
system are: the new system is an improvement over the old, and high level 
vectors offer information beyond that available in the forecast; mid- level 
vectors are of questionable utility and target selection/quality control 
procedures need to be improved; the "slow bias” problem persists and is not 
relieved by the current automatic quality control; finally, manual editing 
of the data enhances the quality. 

NESDIS PROCESSING OF METEOSAT 

Due to the aging of GOES -7 and the launch delay of the next generation 
GOES- I, NESDIS has established a requirement for the processing of CMW from 
the METEOSAT- 3 which has been moved to 50 W /5/. For this processing (and 
also for GOES - 1 ) the 13.7 micrometer measurement used in C02 cloud height 
assignment is not available. Consequently, CIMSS is adapting a variant of 
the height assignment method used at EUMETSAT wherein the infrared window 
temperature representative of opaque cloud is obtained by utilizing 
simultaneous measurements of the 11 micrometer window and the 6.7 micrometer 
water vapor channels /6 /. The success of the method is predicated on the 
fact that the radiances for two spectral bands vary linearly with cloud 
amount. Thus an x/y plot of water vapor radiances versus window radiances 
in a field of varying cloud amount will be nearly linear. This relationship 
is used in conjunction with calculations of outgoing radiance for both 
spectral channels made from the temperature and moisture profiles provided by 
the NMC forecast. A number of calculations are made assuming opaque cloud 
at different levels. These are represented by a curve on the x/y plot, and 
the intersections of curve and line occur at clear sky and opaque cloud 
radiances. The latter, converted to a radiance temperature , is matched to 
the forecast temperature profile to estimate cloud pressure. Fig. 2 shows an 
application of this process. Due to measurement noise and variations in 
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cloud height and emissivity over the 15 X 15 fields of view (centered on the 
cloud target) the radiance pairs do not form a straight line, and in the 
CIMSS application a clustering algorithm is used to define "warm" and "cold" 
cluster centroids which then define the line. 

Pig 2 

Preliminary statistical comparisons of the intercept and C02 methods 
done on the same sample of GOES -7 measurements have shown discrepancies, but 
these are close to the expected inaccuracy of the C02 method estimated to be 
about 50 HPa. The intercept method has also been compared with the infrared 
window estimate (used in the earlier manual processing system) and is clearly 
closer to the C02 assignment /7/. Comparisons of CMW using the intercept 
method have also been done with collocated rawinsondes with the results shown 
in Table 4. The accuracy is slightly worse than that achieved with the GOES- 
7, though the sample size is minimal for drawing any firm conclusion. Also, 
the bias to the forecast (CMV-G) is accentuated, because the height 
assignment is more frequently changed as will be discussed below. This 
problem is probably best alleviated by increasing the density of the CMW 
targets which in turn will lead to more winds for better consistency 
checking, reducing the dependence on the forecast. It would be desirable to 
see magnitudes which are similar between CMV-R and CMV-G. Our general 
conclusion is that CMW are slightly degraded with the substitution of the 
intercept for the C02 method, even with the inclusion of the 
assimilation/quality control. 

WATER VAPOR WINDS 

Winds derived from water vapor imagery have been generated experimentally for 
a number of years at CIMSS/NESDIS . Motivation for this pursuit has been the 
belief that motions can be obtained where no cloud is visible, and that these 
motions represent mid- level flow not well represented in infrared window 
measurements. Both of these beliefs have been controversial. Motions can be 
derived in "clear" air, but these have been demons trably less accurate than 
the infrared window vectors, and EUMETSAT has determined that vectors from 
the water vapor imagery are only reliable when clouds are in the target. 
Under this constraint, the resulting vectors are generally not mid- level /4/. 
However, recent work at CIMSS with the METEOSAT 3 water vapor channel (with a 
better horizontal resolution than the GOES) somewhat revitalizes hope. Fig. 3 
shows an example of coverage obtained for a recent case when the METEOSAT 
water vapor channel was used at CIMSS to provide winds to the National 
Hurricane Center during the approach of hurricane Andrew. Many examples are 
apparent where winds have been obtained (and passed quality control) in clear 
air. Most of these, as discussed below, were assigned to mid levels, and 
collocated rawinsondes for a sample of 39 on this date gave vector rms errors 
of 6.17 and 6.85 ms-1 for the water vapor winds and the NMC forecast 
respectively (with a mean rawinsonde speed of 12.1 ms-1). These statistics 
are competitive with those shown in Table 1 for the GOES infrared window 
winds . 

Fig 3 

It is instructive at this point to consider the height assignment applied to 
the three wind types discussed above: C02 for the GOES infrared window; 
intercept for the METEOSAT infrared window; and infrared histogram for the 
METEOSAT water vapor. Each is subject to revision by the 
assi m ilation/quality control step, and the result for (one day's processing 
of) each is shown in Fig. 4. The C02 heights of the infrared window vectors 
are only occasionally changed in the assimilation step, and a correlation of 
.98 exists between original and final assignment. The intercept method 
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heights are changed somewhat more, leaving a correlation of .80. The 
histogram heights are changed the most, leaving a correlation of only .64. 
These results support the anticipated accuracies of the three height 
assignment methods, and corroborate those reported by Nieman et al. I’ll . It 
is of interest to note that the adjustments for the intercept and histogram 
methods considerably increase the variance of the assignment, with the 
standard deviation increasing from 73 to 111 HPa and from 39 to 60 HPa 
respectively. The smaller variance of the water vapor vector final 
assignment of heights is again anticipated, and it is encouraging to see a 
large number of vectors, initially assigned near 400 HPa, brought down to 
mid- levels . 

Pig 4 


GRADIENT WIND ESTIMATES FROM VAS 

The gradient winds which are provided to the NMC are produced at 4 levels 
(850, 700, 500 and 400 HPa) over an area bounded by 23 • 49 N and 40 -160 W. 
Winds are attempted at a density of every 20 VAS fov (but using the standard 
11 x ll fov averaging) or on a grid mesh of approximately 160-170 km. An 
example of coverage is given in Fig. 5. Success is naturally predicated on 
the absence of cloud. Also, normal quality control is performed on the VAS 
tesperature retrievals, and failing retrievals do not participate in further 
processing. Successful retrievals are hydrostatically integrated from an 
objective analysis of the 1000 HPa height surface to produce geopotential 
heights at the standard pressure surfaces. These are then processed using 
the 2 -dimensional version of the recursive filter objective analysis model 
/8/ to provide geopotential thickness analyses for 1000-700, 700-300, and 
300-100 EPa layers. The 12 -hour NMC forecasts (Aviation) of 
temperature/dewpoint which are used to provide the first guess profiles for 
the retrievals, are similarly integrated on a regular grid to provide 
background fields for the analyses. The fit of the VAS data to these 
thickness analyses provides a second level of quality control, and 
measurements which deviate by 15, 25, and 20 meters, for the three layers 
respectively, are flagged as erroneous. The remaining reports, and the 
background fields, are reanalyzed to provide geopotential surfaces at the 
850, 700, 500 and 400 HPa levels. Analyses are produced on a 1.2 degree 
pseudo Mercator mesh and are quite heavily smoothed. 

Gradient winds are calculated using derivatives obtained from a least squares 
fit of a polynomial geopotential surface to the 25 points centered on the 
grid point closest to each successful VAS retrieval. (Note, the gradient is 
not representative of the actual retrieval location but may be as much as .8 
degrees away.) This is also a smoothing step, in the sense that the 
gradient represents a surface over about 6*6 degrees. The gradient winds 
are finally examined by a meteorologist who further edits before transfer to 
the NMC. 

The gradient winds are routinely evaluated in collocation statistics at 
NBSDIS. Results obtained during August 1992 are given in Table 5. 

This table shows that the accuracy of the gradient wind is comparable to that 
of the CMW, and again, the satellite product is able to add information to 
the forecast fields. In this case we can further show that the satellite 
apparently adds information to an analysis which includes the height reports 
of the same rawinsondes used for verification. The reason for this is 
probably that the satellite data are more densely distributed. 
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TABLE 5 

Experiments have been conducted at CIMSS and NMC introducing the VAS gradient 
winds, as observed winds, into data assimilation models. The effect has been 
negligible. Modest changes to the wind fields at large scales are quickly 
lost. Also, there are occasions when the gradient wind is a poor 
approximation to the actual flow and the data are detrimental to a wind 
analysis. An example can be seen in the upper mid-west in Fig. 5. One might 
then ask why the VAS data should not be introduced as temperatures , allowing 
the assimilation to extract the gradient. The difficulty with such an 
approach is that we have found that rather heavy smoothing of the 
geopotential analysis is required if the VAS gradient winds are to be as 
effective as shown in Table 5. The smoothing is necessary because the data 
are noisy. We have not, to date, been able to show, consistently, additional 
information over that contained in the NMC forecasts when temperature (or 
geopotential) is compared to radiosondes. It is therefore not obvious that 
the VAS can be beneficially incorporated, as temperatures , with other 
temperature data. In view of this, it is probably optimum to use the VAS 
data as gradient winds to shape the geopotential analyses, but not to use the 
gradient winds as observed winds because of the ageostrophic regimes. We 
shall be pursuing this approach at CIMSS. 

SUMMARY 

Four types of winds generated from satellite data have been considered. The 
first of these, the new operational winds from GOES infrared window 
measurements, is shown to be a notable improvement over the previous manual 
wind generation. High level vectors are shown to contain information beyond 
that available in the NMC 12 hour forecast. However, a long standing speed 
bias remains. Infrared winds from the METEOSAT are under development, 
utilizing a two channel "intercept" height assignment similar to that 
employed at EUMETS AT . While this assignment is clearly superior to older, 
histogram methods, the winds appear to be slightly inferior to GOES processed 
with an the C02 height assignment, even after the assimilation/quality 
control step. Water vapor winds generated from the METEOSAT are regenerating 
hope for a source of mid- level motions. Again these rely on the 
assimilation/quality control. Finally, gradient winds derived from VAS 
temperature profiles have been evaluated, and shown to contain information 
over the forecast at 850, 700, 500 and 400 HPa. These should complement the 
higher level CMW in cloud free areas. 
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Table l. February- -July, 92 statistics of automated CMW versus rawinsondes 
(R) and NMC 12 hour forecast guess (G) . hi refers to <400 HPa. mid refers to 
699<P<399 HPa. Units are in ms-l. 


RMS Vector Error 
(sample) 

G-R 

CMV-R 

CMV-G 

hi (1163) 

7.8 

7.5 

5.8 

mid (337) 

5.9 

6.3 

4.0 

Mean Speed 

G 

CMV 

R 

hi 

20.2 

21.7 

23.0 

mid 

15.9 

12.7 

14.3 


Autoeditor Only (Ho Manual) 


RMS Vector Error 

G-R 

CMV-R 

CMV-G 

(sample) 

hi (1458) 

7.8 

7.8 

5.9 

mid (552) 

5.5 

6.0 

4.0 

Mean Speed 

G 

CMV 

R 

(m/s) 

hi 

18.9 

19.4 

21.2 

mid 

15.2 

11.6 

13.1 


Table 2. February- -July, 91 statistics of manual CMW versus rawinsondes (R) 


RMS Vector Error 

G-R 

CMV-R 


(sample) 

hi (1024) 

9.2 

8.6 


mid (228) 

6.5 

6.3 


Mean Speed 

G 

CMV 

R 

(m/s) 

hi 

12.9 

16.6 

17.2 

mid 

8.1 

10.0 

10.9 
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Table 3. Speed bias, as a function of rawinsonde wind speed for 12-hr 
aviation forecast (AVN) and the processed forecast used in CMSf generation 
(CMWF) . Units are ms-1. 


Speed 

AVN 

CMWF 

0-10 

1.8 

-0.2 

10-20 

-0.0 

-1.8 

20-30 

-0.4 

-2.7 

30-40 

-2.5 

-4.9 

40-50 

-3.5 

-5.6 

>50 

-9.2 

-12.5 


Table 4. Accuracy estimates of CMW processed from METEOSAT - 3 infrared window 
imagery. Sample is for 26-31 August, 1992. 


RMS Vector Error 

G-R 

CMV-R 


CMV-G 

(sample) 






hi (93) 

8.4 

9.4 


5.2 


mid (33) 

6.7 

7.5 


3.8 


Mean Speed 

G 

CMV 


R 


hi 

22.5 

20.8 


25.0 


mid 

12.2 

12.1 


13.0 


Table 5. August 92 statistics of NESDIS VAS gradient winds versus 

rawinsondes (R) . Statistics calculated 

for gradient winds 

generated 

12 hour forecasts (G) and 

radiosonde analyses 

(A) are 

also given. 

RMS Vector Error 

G-R 

VAS-R 


A-R 


850 

5.5 

4.8 


5.1 


700 

4.6 

4.2 


4.3 


500 

5.4 

4.9 


5.0 


400 

6.8 

6.0 


6.4 


Mean Speed 

G 

VAS 

A 


R 

(m/s) 






850 

4.6 

4.7 

4.4 


6.1 

700 

5.7 

5.4 

5.4 


5.9 

500 

8.7 

8.5 

8.3 


8.3 

400 

10.9 

10.3 

10.- 

4 

10.7 
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Fig. 1. CMW bias (CMW-Rawinsonde) as a function of rawinsonde speed. Crosses 
and dashed regression line represent CMW, squares and solid line represent 
NMC aviation forecast. 
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Fig. 2 Example of the "intercept" method for obtaining the temperature of 
opaque cloud from multiple partly cloudy views. See text for explanation. 
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Fig* 3. Coverage of water vapor winds obtained in the vicinity of Hurricane 
Andrew using the imagery from METEOSAT-3 for 23 August, 1992. 




Fig. 4. Comparison of pressure heights resulting from assimilation/quality 
control with three methods used for GOES -7 infrared window CMW (C02), 
METEOSAT-3 infrared window CMW (Intercept), and METEOSAT-3 water vapor CMW 
(histogram) . 
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Fig. 5. VA S gradient wind coverage for 25 August 1992. Example shown is for 
400 HPa. Bold vectors are rawinsonde reports. 
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Abstract 

This paper describes the results from a collaborative study 
between the European Space Operations Center, the European 
Organization for the Exploitation of Meteorological Satellites, the 
National Oceanic and Atmospheric Administration, and the Coop- 
erative Institute for Meteorological Satellite Studies investigating the 
relationship between satellite-derived monthly mean fields of wind 
and humidity in the upper troposphere for March 1994. Three 
geostationary meteorological satellites GOES-7, Meteosat-3, and 
Meteosat-5 are used to cover an area from roughly 1 60°W to 50°E. 
The wind fields are derived from tracking features in successive 
images of upper-tropospheric water vapor (WV) as depicted in the 
6.5-p absorption band. The upper-tropospheric relative humidity 
(UTH) is inferred from measured water vapor radiances with a 
physical retrieval scheme based on radiative forward calculations. 

Quantitative information on large-scale circulation patterns in the 
upper troposphere is possible with the dense spatial coverage of the 
WV wind vectors. The monthly mean wind field is used to estimate 
the large-scale divergence; values range between about -5 x 10 * 
and 5 x 10* sec 1 when averaged over a scale length of about 
1000-2000 km. The spatial patterns of the UTH field and the 
divergence of the wind field closely resemble one another, suggest- 
ing that UTH patterns are principally determined by the large-scale 
circulation. 

Since the upper-tropospheric humidity absorbs upwelling radia- 
tion from lower-tropospheric levels and therefore contributes signifi- 
cantly to the atmospheric greenhouse effect, this work implies that 
studies on the climate relevance of water vapor should include three- 
dimensional modeling of the atmospheric dynamics. The fields of 
UTH and WV winds are useful parameters for a climate-monitoring 
system based on satellite data. The results from this 1 -month 
analysis suggest the desirability of further GOES and Meteosat 
studies to characterize the changes in the upper-tropospheric mois- 
ture sources and sinks over the past decade. 


‘European Organization for the Exploitation of Meteorological 
Satellites, Darmstadt, Germany. 

*NOAA/NESDIS, Madison, Wisconsin. 

“Cooperative Institute for Meteorological Satellite Studies, University 
of Wisconsin — Madison, Madison, Wisconsin. 

° European Space Agency, European Space Centre, Darmstadt, 
Germany. 

Corresponding author address: Dr. Paul Menzel, NOAA/NESDIS, 
1225 West Dayton St., Madison, Wl 53706. 

In final form 26 January 1995. 

©1 995 American Meteorological Society 


1. Introduction 

The water vapor band of the geostationary satel- 
lites is routinely used for extracting upper-tropospheric 
humidity (UTH) and water vapor (WV) wind fields. 
Apart from the utility of both products in numerical 
weather prediction and model validation, it has be- 
come evident that UTH and WV winds constitute an 
important source of climatological information. Only a 
geostationary platform provides measurements of 
water vapor amount and motion on a large scale from 
a single observing system; thus, these measurements 
serve as truly independent datasets for comparisons 
with other climatologies. Through joint efforts at the 
European Space Operations Center (ESOC) and the 
Cooperative Institute for Meteorological Satellite Stud- 
ies (CIMSS), data from three geostationary satellites 
(GOES-7 at 112°W, Meteosat-3 at 75°W, and 
Meteosat-5 at the Greenwich Meridian) are processed 
for March 1 994 with similar algorithms to derive UTH 
and WV winds covering 160°W-50°E (Fig. 1). Monthly 
mean observations are computed from two UTH and 
WV wind fields per day (1200 and 0000 UTC) for 
GOES-7 and Meteosat-3 ; for Meteosat-5 eight UTH 
and four WV wind fields are used per day. 


2. Importance of UTH 

Satellite determinations of UTH have been shown 
to be useful for deriving climatologies (van de Berg et 
al. 1991; Wu et al. 1993) that can be compared with 
those derived from models (Del Genio et al. 1994; 
Stephens et al. 1993; Soden and Bretherton 1993; 
Schmetz and van de Berg 1994). The comparisons 
reveal significant differences between UTH from sat- 
ellite observations and those in general circulation 
models. 

Figure 2 illustrates that all altitude levels of the 
watervaporprofile, including those represented by the 
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3. UTH algorithms 



Fig. 1 . Coverage from the three geostationary satellites 
{GOES-7 at 1 12°W, Meteosat-3 at 75°W, and Meteosat-5 at the 
Greenwich Meridian) out to 50° from nadir. 


UTH, are important to the water vapor greenhouse 
effect. The greenhouse effect (e.g., Inamdar and 
Ramanathan 1994) is defined here as the difference 
between the longwave flux emitted by the surface and 
the outgoing radiative flux (OLR) escaping from the 
top of the atmosphere. A radiative transfer model is 
used to compute the change of the OLR as a function 
of perturbations in the vertical water vapor profile. A 
10% increase in relative humidity is added to a stan- 
dard tropical profile at a single level at a time; the levels 
are at 1-, 3-, 5-, 7-, 9-, 11-, and 13-km altitudes. The 
sensitivity is expressed as the decrease in OLR di- 
vided by the increase in total precipitable water, 
del(OLR)/del(PW). Figure 2 shows that this sensitivity 
increases from about 2.8 to about 1200 W m 2 (g cm 2 ) 1 
as the 10% relative humidity perturbation is shifted 
from 1- to 13-km altitude (solid curve). A mixing ratio 
profile typical for a tropical atmosphere is also shown 
in Fig. 2 (dashed curve). The greenhouse effective- 
ness of a single H 2 0 molecule (solid curve) rapidly 
increases with altitude. This is to be expected since 
radiation from the lower atmosphere absorbed by H 2 0 
molecules high in the troposphere is reemitted at 
colder temperatures, in contrast to reemission at warm 
temperatures by the lower-atmospheric H z O molecules. 
The exponential decrease of the water vapor with 
altitude balances the increase of the greenhouse 
effectiveness of a single water vapor molecule. The 
combined effect suggests that the sensitivity of OLR to 
perturbations in the water vapor profile in a tropical 
atmosphere is roughly constant up to an altitude of 
about 10 km. Thus, Fig. 2 indicates the need for 
observations of the UTH, as well as observations of 
total precipitable water (which is dominated by low- 
level water vapor) in studies of the atmospheric green- 
house effect. 


a. Meteosat-5 

The UTH is inferred from the W V band (5. 7-6. 9 n), 
which has a spatial resolution of 5 km x 5 km. The WV 
band observes the earth’s atmosphere in a spectral 
region characterized by strong H 2 0 absorption and is 
sensitive to approximately the upper 3 mm of total 
water vapor (Schmetz and van de Berg 1994). In a 
cloud-free atmosphere, this WV band senses typically 
down to a pressure level of 500-600 hPa. The height 
allocation of the broad contribution function depends 
on the actual profiles of temperature and humidity and 
the satellite viewing angle (Fischer et al. 1981). 

The UTH retrieval follows the approach described 
by Schmetz and Turpeinen (1988). Radiative forward 
calculations are performed for different relative hu- 
midity profiles. The necessary temperature profiles 
are taken from short-term model forecasts provided 
by the European Centre for Medium-Range Weather 
Forecasts. The different relative humidity profiles are 
constructed as follows: below 600 hPa the forecast 
profile is used, between 600 and 300 hPa the relative 
humidity is fixed at one of six values (1%, 5%, 10%, 
20%, 40%, 100%), and above 300 hPa the relative 
humidity decreases linearly to 0% at 100 hPa. Com- 
parisons of the UTH retrieval (van de Berg et al. 1 991 ) 
with radiosondes show a root-mean-square differ- 
ence of less than 10%. The Meteosat-5 UTH calcula- 
tions use clear-sky radiances obtained from a multi- 
spectral image classification for segment areas 
(32 x 32 pixels), ranging from 160 up to 230 km on a 
side. 


dOLR / dPW [ W m ’/(jc m - ’)] Mixing Ratio (g/kg) 



Altitude (km) 


Fig. 2. Sensitivity of the outgoing longwave radiative flux density 
(OLR in W nr 2 ) to changes in precipitable water (PW in g crrr 2 ). PW 
has been varied at different levels, and the sensitivity is expressed 
as the decrease in OLR over the increase in PW. The dashed curve 
corresponds to the mixing ratio profile characteristic of a tropical 
standard atmosphere. 
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The Meteosat-5 UTH in this study is somewhat 
different from that presented in van de Berg et al. 
(1991). That work confined valid UTH retrievals to 
segment areas that were totally free of medium- and 
high-level cloud so that the UTH was representative of 
the whole segment. Here the UTH retrieval is extended 
to segments that are partially covered with high- and 
medium-level cloud by using the WV radiances from 
the portion of the segment that contains clear sky. This 
tends to moisten the UTH, especially in the tropical 
convective regions where high-level clouds prevail. 

The operational calibration of the Meteosat WV 
band is based on radiative transfer calculations with 
radiosonde profiles and collocated clear-sky WV ob- 
servations. A recent recalibration algorithm has been 
added for this study, which improves calibration per- 
formance (Schmetz and van de Berg 1994). 

b. Meteosat-3 

Meteosat-3 data are similar to those of Meteosat-5, 
differing only in spectral range (5. 8-7. 3 p). However, 
the UTH retrieval from Meteosat-3 differs from the 
Meteosat-5 UTH retrieval in several aspects. The 
major difference is that the Meteosat-3 UTH is re- 
trieved for each pixel regardless of cloud conditions, 
so it has much higher spatial resolution . When the WV 
radiance for a given pixel is affected by cloud, the UTH 
does not represent a column-mean relative humidity 
and can indicate up to 100% relative humidity. Includ- 
ing cloudy pixels in the retrieval increases the monthly 
mean UTH by about 10% (absolute) in convective 
areas; no change is observed in areas without high- 
level cloud. As a result, monthly mean Meteosat-3 
UTH values are slightly more moist than those from 
Meteosat-5 in cloudy regions. 

There are also several minor algorithm changes. 
Soden and Bretherton (1993) showed that, under 
normal conditions, the WV band brightness tempera- 
ture is close to a logarithmic function of the UTH. 
Accordingly, the forward calculation of UTH from 
Meteosat-3 , which is used to generate the lookuptable 
for WV vapor band brightness temperature, is now 
performed on a different set of values (1%, 5%, 15%, 
30%, 60%, 100%) and the interpolation is changed 
from linear to loglinear. Finally, the CIMSS forward 
model is used for radiative transfer computation. 

c. GOES-7 

The GOES-7 WV band (6.4-7. 1/j) has a spatial 
resolution of 7 km x 7 km at nadir. While the UTH 
retrieval algorithm used with GOES-7 is identical to 
that used with Meteosat-3, the UTH fields from the 
two satellites are somewhat different because the two 
instruments sense in a slightly different part of the 
spectrum. In a clear atmosphere, the GOES-7 WV 


band typically senses down to a pressure level of 
400-500 hPa, depending on view angle and atmo- 
spheric moisture content (Hayden 1 988). Thus, com- 
pared to the Meteosats, the GOES-7 WV band senses I 
radiation from a slightly higher layer in the upper , 
troposphere. For March 1994 over an area (30°N- 
30°S, 82°W-105°W) with relatively equal and small ' 
viewing angles for both GOES-7 and Meteosat-3, 
tests showed the GOES-7 UTH was typically 5% 
higher than the Meteosat-3 UTH. 

The empirical distribution function (EDF), described 
by Weinreb et al. (1989), is employed to eliminate this 
intersatellite bias. For each of the concurrent observa- 
tions from the above-mentioned area (62 for March), 
an EDF is constructed that is used to adjust the 
GOES-7 UTH to have the same statistical distribution 
as the Meteosat-3 UTH. The unadjusted GOES-7 
UTH, which is representative of a higher layer in the 
uppertroposphere, is consistently more moistthan the 
Meteosat-3 UTH. This suggests that the relative hu- 
midity for the period and area under consideration 
increases with height in the upper troposphere. Al- 
though such profiles are not uncommon, it is generally 
agreed that relative humidity decreases with height for 
longer-term averages. A possible explanation is that 
Meteosat-3 calibration was found to be warm with 
respect to GOES-7, which would make the UTH drier 
(Menzel et al. 1 994). Further investigation is required 
to clarify this issue. 

d. Mean UTH for March 1994 

Figure 3 shows the monthly mean UTH for March 
1994 from all three satellites. The adjusted GOES-7 
UTH merges very well with the Meteosat-3 UTH with- 
out noticeable discontinuity, while the transition from 
Meteosat-3 to Meteosat-5 is distinguishable (prima- 
rily because different UTH retrieval procedures are 
used). A unified retrieval for all satellites would likely 
alleviate these boundaries. 

Highest relative humidities of about 60%-70% are 
observed over the tropical convective zones of South 
Africa and South America. The relative humidity val- 
ues are computed with respect to water vapor satura- 
tion pressure. (At an altitude of 400 hPa, the ratio of 
saturation pressure over ice to that over water is about 
0 . 8 .) 


4. WV wind fields 

a. Meteosat-5 

The method for deriving Meteosat-5 WV winds at 
ESOC is described by Laurent (1993) and Holmlund 
(1993). The height assignment of the WV winds fol- 
lows two approaches, (a) If the target area for tracking 
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Fig. 3. The monthly mean UTH for March 1994 from all three satellites. Values have been averaged to roughly a 1° lat-long grid. The 
adjusted GOES-7 UTH merges very well with the Meteosat-3 UTH without noticeable discontinuity. On the other hand, the merge line 
between Meteosat-3 and Meteosat-5 \s distinguishable, primarily because different UTH retrieval procedures are used. 


the displacement contains cloudy features, the height 
assignment is the same as for the high-level infrared 
cloud motion winds (Schmetz et al. 1993). Opaque 
clouds are assigned to the pressure level where the 
brightness temperature fits a short-term forecast pro- 
file; semitransparent clouds are assigned to a pres- 
sure level using the IR and WV bands, (b) If the target 
area is cloud free, the coldest 25% of the WV pixels in 
a target area is used to compute a brightness tempera- 
ture for the height assignment (Laurent 1993). The 
brightness temperatures are converted to pressure 
heights via short-term forecast temperature profiles. 
The adequacy of this approach has been verified in 
detailed comparisons with rawinsondes (Holmlund 
1993). 

b. Meteosat-3 and GOES-7 

Water vapor wind sets derived at CIMSS from 
Meteosat-3 and GOES-7 are processed from hourly 
full-disk images using fully automated procedures 
[see Merrill et al. (1991) for further details]. Potential 
targets are identified by maximum brightness tem- 
perature gradients exceeding predetermined thresh- 
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olds within a selected target region. Heights are as- 
signed using the WV band brightness temperature of 
the coldest 25% of the target region matched to a 
short-term forecast profile. The vertical distribution 
of the pressure heights is typically in the range of 
200-500 hPa with a peak near 300 hPa. As a final 
step, an objective quality control procedure is invoked 
once the wind vector field has been derived (Hayden 
and Velden 1991). Statistical comparisons between 
WV wind vectors derived at CIMSS and collocated 
Western Hemisphere rawinsondes show rms differ- 
ences that are close to those of current operational 
upper-air cloud-drift winds (Velden et al. 1 994). Water 
vapor winds have also been shown to have positive 
impact on NWP (Velden 1995). 

c. Mean WV wind field for March 1994 

The daily WV winds produced by both CIMSS 
(2 d - ') and ESOC (4 d ’) for March 1994 are used to 
calculate a monthly mean field by averaging all wind 
vectors contained within 2.0° lat-long geographical 
regions (grid points). A monthly mean gridpoint motion 
vector value is considered to be valid when the num- 
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Fig. 4. The mean WV wind vectors displayed in kts for March 1 994. (Vectors every 4° are shown.) Blue indicates vectors with pressures 
greater than 400 hPa; brown indicates vectors with pressures between 300 and 399 hPa; and green corresponds to vectors with pressures 
less than 300 hPa. 


ber of observations exceeds 20. (Only grid points near 
the limb fail this criterion.) The monthly mean gridpoint 
height value is calculated by averaging individual 
vector height assignments. Figure 4 illustrates that the 
mean WV wind vectors show good consistency and 
the heights are confined to adequately describe the 
upper-tropospheric layer mean flow. 


5. UTH and WV wind field divergence 

To associate the observed upper-tropospheric cir- 
culation with the observed UTH, divergence fields are 
computed with a finite-difference scheme from 
the monthly mean WV winds. The dashed contours in 
Fig. 5 indicate mean divergent areas that are often 
associated with water vapor emerging from source 
regions (deep convective areas). For example, the 
pattern in the UTH over South America is clearly 
correlated with the divergence field. 


The solid contours in Fig. 5 correspond to upper- 
tropospheric convergence, most prominent over sub- 
tropical subsidence regions. The observed convergence 
values appear quite reasonable: assuming that the 
convergence takes place in a layer of about 1 50 hPa, 
the associated subsidence is about 60 hPa d _1 , which 
compares favorably to the observed subsidence in the 
trade wind regions (Betts and Ridgway 1989). 

A comment is warranted concerning the height 
variation of the WV winds. Inspection of Fig. 4 shows 
that the monthly mean altitudes of the WV winds range 
approximately between 250 and 350 hPa. The spatial 
changes are smooth, and the standard deviation at a 
geographical location is only 30-70 hPa. Thus, the 
computed divergence/convergence patterns do not 
correspond to a constant pressure surface but to a 
smoothly varying surface with higher altitudes in the 
convective area than in the subsidence regions. 

One should not expect a perfect match between 
UTH and wind field divergence. Over tropical areas 
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Fig. 5. Simultaneous display of the upper-level convergence fields (derived from WV winds) and the observed UTH for March 1 994. The 
dashed contours correspond to divergent areas in the upper troposphere; solid contours represent convergence. Contours are shown in 
intervals of 2 x 10 6 s ’. The UTH is displayed as a color image with a relative humidity scale displayed. (Blue is dry, green is in the middle, 
and red is moist.) 


especially, the UTH maxima are associated with con- 
vection usually forced by surface heating and/or con- 
vergence. On the other hand, the WV wind field is 
primarily representative of upper-level large-scale flow, 
with secondary influence from large convective re- 
gimes. The UTH and divergence maxima should be in 
close proximity, but they need not be aligned. An 
example is evident in Fig. 5 over Africa, where the UTH 
maximum is shifted to the south of the wind divergence 
maximum. The position of the divergence maximum 
appears to be influenced by the entrance region of the 
subtropical jet over Africa. However, over the South 
American continent, where no jet stream is observed, 
the divergence maximum matches well with the UTH 
maximum. 

The congruence of the UTH and WV wind diver- 
gence patterns suggests a physical relationship in 
which the UTH patterns are largely influenced by 
upper-tropospheric transport. Picon and Desbois 
(1990) indicate that the upper-tropospheric humidity 
field is chiefly determined by large-scale dynamics. 
The clear-sky UTH climatology of van de Berg et al. 
(1991) corroborates this statement, as the geographi- 
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cal distribution of monthly UTH fields is reminiscent of 
the Hadley and Walker circulation patterns. Soden 
and Bretherton (1993) also point out that the lower 
troposphere is the prime source of upper-tropospheric 
water vapor and that deep convection is the transport 
mechanism. However, the novel message of this 
study is that both UTH and the upper-tropospheric 
dynamics can be inferred simultaneously from geo- 
stationary satellites. As such, the satellite-derived 
UTH and WV winds constitute a valuable source for 
climatological diagnostics and for comparison with 
large-scale circulation models. 


6. Summary 

Monthly mean fields of the upper-tropospheric rela- 
tive humidity (UTH) and water vapor (WV) winds from 
Meteosat-3, Meteosat-5, and GOES-7 WV band ob- 
servations have been analyzed for March 1 994. The 
fields provide a concise description of monthly mean 
upper-tropospheric moisture and motion fields. The 
scope of the present study is a qualitative demonstra- 
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tion that UTH and upper-tropospheric dynamics are 
directly observable on a nearly global scale; a future 
quantitative analysis will include a much larger dataset. 

The comparison of the March 1994 UTH and WV 
wind fields shows the influence of atmospheric dy- 
namics on the distribution of the upper-tropospheric 
moisture. Since upper-tropospheric moisture is an 
important factor in the atmospheric greenhouse ef- 
fect, the satellite-observed fields of UTH and WV wind 
are important for climate model validation. Most of the 
previous comparisons of climate models with ob- 
served moisture fields used only total precipitable 
water vapor observations over the oceans. Since low- 
level moisture dominates precipitable water, such 
comparisons do not validate the upper-tropospheric 
moisture fields in models. 

This study using satellite measurements docu- 
ments an example where the upper-tropospheric 
moisture distribution is strongly influenced by the 
large-scale dynamics. Longer-term datasets of both 
UTH and WV winds covering seasonal and interannual 
variations would be useful for comparison with climate 
models and for validation of the dynamical aspects of 
the models. Such a study could help reveal how the 
large-scale circulation would be affected by climate 
change. 

This study also suggests that use of WV winds 
should improve representation of upper-tropospheric 
moisture fields in numerical models. This may be 
important to medium-range forecasts since the radia- 
tive heating/cooling critically depends on upper-tropo- 
spheric moisture. The importance is amplified by the 
fact that the moisture field influences the cloud field, 
which has an even larger effect on radiative heating/ 
cooling and also implies energy transformations 
through phase changes. 

Both UTH and WV winds could be very desirable 
contributors to a climate-monitoring system. The moni- 
toring of both quantities will soon be possible on a 
nearly global basis since future geostationary satel- 
lites from China and Japan will also be capable of 
radiance measurements in the WV band. Further 
detailed studies are needed to fully understand the 
features in simultaneous UTH and WV wind fields. A 
continuous operational derivation of the two products 
would provide the basis for investigation of seasonal 
and interannual variations. 
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Summary 


Tropospheric motions deduced from sequential water vapor imagery provided by 
geostationary meteorological satellites can be utilized to infer wind fields. Temporal tracking of 
moisture features yields spatially coherent vector fields in both cloudy and cloud-free regimes. 
These observations can be employed to augment existing operationally-available data in order 
to provide improved upper-tropospheric wind analyses over meteorological scales ranging from 
sub-synoptic to global. It is demonstrated through assimilation of these data into numerical 
weather prediction systems that modest forecast improvements can be realized. Forecast 
impact experiments yield small but positive results in: 1) hurricane track forecasts, 2) regional- 
scale prediction, and 3) global-scale prediction. Ongoing research and future prospects are 
discussed. 
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1. INTRODUCTION 


Water vapor imagery from geostationary meteorological satellites has been available for 
over a decade. These data are employed extensively by operational analysts and forecasters, 
mainly in a qualitative mode (Weldon and Holmes 1991). Early attempts to quantify the data by 
employing manual methods to track moisture features in animated image sequences met with 
modest success (Stewart et al. 1985; Hayden and Stewart 1987). More recent studies have 
demonstrated the feasibility of automated techniques to produce wind vector fields comparable 
in quality to operational cloud-tracked winds (Laurent 1993; Velden 1993; Holmlund 1993). 

Development of automated wind-extraction algorithms at the Cooperative Institute for 
Meteorological Satellite Studies (CIMSS) has progressed to the stage of quasi-operational 
testing and implementation (Nieman et al. 1995). In collaboration with the NOAA National 
Environmental Satellite Data and Information Service (NESDIS), efforts at CIMSS have 
focused on the research, development and demonstration of these algorithms towards the goal 
of acceptance by the National Meteorological Center (NMC) for the routine inclusion of the 
data into the operational data stream. In this paper, results pertaining to the assimilation of the 
water vapor tracked winds into selected numerical weather prediction (NWP) systems is 
presented. The impact on varying meteorological regimes is evaluated by examining the 
influence of the data on numerical prognoses ranging from hurricane to global scales. 

2. DATA SET DESCRIPTION 

Wind vectors derived from water vapor image sequences are generated by automated 
algorithms following the same general philosophy as with cloud tracking (Merrill et al. 1991). 
Potential targets (moisture features or structure) are first identified in an image. Employing 
correlation techniques, an attempt is then made to track the targets in subsequent images at 
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hourly intervals. Successfully tracked targets (vectors) are then assigned heights based on an 
optimal match between target radiance (brightness temperature) and a guess temperature 
profile (usually a global forecast). In most cases, this assigned height is representative of a 
layer of moisture which is being tracked, since the pixel radiances are achieved through 
contributions of energy over the atmospheric depth represented by the attendant weighting 
function. It has been clearly demonstrated for future considerations that full disk (global view) 
wind sets can be created using these algorithms on a time scale commensurate with real time 
operations. For further details of this entirely automated procedure see Merrill et al. 1991, and 
Nieman et al. 1995. 

As a final step, an objective quality control (QC) procedure is automatically invoked once 
the wind vector field has been completed (Hayden and Velden 1991). This procedure involves 
three principal steps: 1) A preliminary objective 3-dimensional analysis of the wind data which 
includes any ancillary data and a guess background field, 2) A re-evaluation of the originally 
assigned vector heights, and 3) An objective re-analysis using the vectors with re-assigned 
heights and the preliminary analysis as the guess background field, which provides final quality 
flags. This entire procedure is shown schematically in Figure 1. The objective analysis 
employed is of the successive approximation type, and is basically an adaptation of the 3- 
dimensional recursive filter described by Hayden and Purser (1988). The analysis is performed 
over a domain dictated by the vector coverage. The NMC global assimilation system 6 hour 
forecasts are used for background fields. 

At each iteration of the preliminary analysis, a quality indicator is assigned to each 
observation based on a simple empirical function which basically represents the fit of the data 
relative to neighboring observations and the background. At this point individual vectors are 
given initial quality flags based on their fit to the analysis. The second step involves 
reconsidering the initial height assignment which was loosely based on the attendant 
weighting function and resulting brightness temperature. For each vector, an attempt is made 
to minimize, in the vertical, the following penalty function: 
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((V(p) - SDWV)/WV)2 + (( T (p) . TC)/WT)2 + ((P - PW)/WP)2 


where SDWV is the satellite-derived wind vector; TC and PW are the originally assigned 
temperature and pressure of the wind vector, respectively; V(p), T(p) and P are obtained from 
the background fields; and WV, WT and WP are empirically determined weights currently set 
at 2, 10 and 100, respectively. The third term is included to prevent large deviations in re- 
assignment, since we are relatively secure in the fact that the wind heights are fundamentally 
limited by the spectral response function to a layer in the upper-troposphere. The height 
reassignment of an individual vector may fail if the minimum exceeds predetermined 
thresholds. In this instance, the vector is flagged as a bad vector and is not considered in the 
final re-analysis. Evaluation of this QC procedure at CIMSS has shown it to increase the 
number of vectors retained, improve vector field RMS errors (vs. collocated rawinsondes), 
while altering mean field heights (pressure) by less than 25 Mb. While somewhat instrument 
(satellite) dependent, most resulting vectors are assigned final pressure-heights in the 200-500 
Mb range, with a peak near 300 Mb. 

The full analysis is repeated in step 3 using the water vapor vectors with reassigned 
heights, and the output of the first analysis as a background field. The output of this procedure 
provides the final quality flags. Individual vectors which do not attain a predetermined quality 
are again flagged as bad vectors at this stage, and not passed on to the user. In the 
experiments outlined below, this was the case. However, in operational practice, the option will 
exist for the user to obtain and re-examine the vectors flagged as bad. Further details on this 
objective editing procedure can be found in Hayden and Velden (1991). 

Currently, water vapor wind vectors can be produced from both GOES and Meteosat 
satellites. An example of combined data coverage from GOES-7 and Meteosat-3 is shown in 
Figure 2. From a purely qualitative standpoint, the horizontal coverage and vector coherence is 
quite good. The very high-density coverage near the center of Fig. 2 is due to overlap in 
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satellite coverage. It is noteworthy that, in general, there is very good agreement between 
measurements from the two different satellites. This has also been quantitatively confirmed in 
terms of vector RMS errors (Velden et al. 1994). 

Evaluation of the water vapor winds can be accomplished through two approaches: 
statistical validation against collocated rawinsondes, and numerical model impact studies. 
Earlier studies have shown the water vapor tracked wind estimates to be comparable in quality 
to upper-tropospheric operational cloud tracked winds (Laurent 1993), with RMS values in the 
range of 6.5 - 8 m/sec, and speed biases of less than 1 m/sec (Velden 1993; Velden et al. 
1994). The remainder of this paper will focus on the impact of the winds on NWP. 

3. IMPACT ON NUMERICAL WEATHER PREDICTION 

While statistical evaluation of new data types is a necessary step to determine error 
characteristics, the ultimate worth of the observations in a quantitative sense is largely 
determined by their impact on numerical analysis and prediction. In this section, the forecast 
impact realized from numerical prediction systems initialized with the water vapor wind data is 
presented. Three different prediction systems are examined in order to assess the impact of 
the data on various meteorological scales ranging from sub-synoptic (hurricane) to global. 

a) Impact on numerical hurricane track forecasts 

In our first investigation, a hurricane track prediction model (VICBAR, DeMaria et al. 1992) 
is employed to test the sensitivity of the water vapor wind data on numerical hurricane track 
forecasts. VICBAR is a nested, spectral barotropic model that has been run in near-real time at 
the NOAA Atlantic Oceanographic and Meteorological Laboratory Hurricane Research Division 
(AOML-HRD) for the past few years. The initial condition for the forecast model is a vertically 
averaged (mass weighted) deep-layer mean wind over the 850-200 Mb depth of the 
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troposphere, with an added axisymmetric bogus tropical cyclone vortex. The barotropic 
forecast model uses the shallow-water equations, with the forecast storm track determined 
from the location of the relative vorticity maximum on the innermost model mesh. This 
forecasting system has been used to evaluate impacts of other data types (Franklin and 
DeMaria 1992; Velden et al. 1992). 

Water vapor wind sets were produced from Meteosat-3 daily at 1200 UT during several 
North Atlantic tropical cyclones in 1992 and 1993. Winds were derived over a domain covering 
0-50 N, 30-100 W. For the model impact evaluation, experimental forecasts were performed 
based on the following criteria: 1) the storm was within the domain of the generated wind set, 
and 2) the storm intensity was tropical storm strength or better (>35 kts). There were 32 
forecast cases, from 10 different storms (including Hurricane Andrew in 1992) that met these 
criteria. 

Table 1 shows the results of forecasts which were initialized with the water vapor wind sets 
produced at CIMSS plus operationally available data, compared to control runs which were 
initialized on operationally available data only. Considering all of the cases together, the results 
show generally small improvements to the forecasts with the inclusion of the water vapor wind 
data. The improvement in mean track forecast error (MFI in Table 1) ranges from 2.2% at 24 h 
to 6.8% at 72 h. After an adjustment is made for the serial correlation between forecast cases 
(Franklin and DeMaria 1992), none of the improvements were found to be significant at the 
95% confidence level due in part to the small sample size. Another measure of forecast impact 
is the frequency of improved forecasts (FIF), which simply shows the percentage of forecasts 
which resulted in some improvement when the water vapor data were included. The FIF at 72 
h indicates that 71% (15 out of 21) of the forecasts were improved. 

Several of the selected tropical cyclone cases were relatively close to the US mainland. It is 
reasonable to assume that in these cases, the VICBAR model was relatively well initialized by 
the nearby conventional (operational) data base, limiting the potential for the satellite data to 
have a positive impact. On the other hand, storms well out to sea should make better 
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candidates for forecast improvement. To test this hypothesis, the sample was stratified to 
include only those cases east of 70 W. The results of these 23 cases are shown in Table 2. 
While still not statistically significant at the 95% confidence limit, the MFI percentages and FIF 
are a notable improvement over the statistics presented in Table 1. 

Using the stratified sample, it is of interest to examine the distribution of forecast 
differences between the experimental and control runs in order to fully appreciate the FIF 
results. Table 3 shows the distribution of the 72 h forecast differences. Of the eighteen cases, 
only two experimental forecasts were notably degraded (defined here as greater than 100 Km) 
relative to the control forecasts. One of these poor forecasts was diagnosed as a case in which 
the storm did not follow the deep layer flow. Rather, it followed a shallow-layer flow, in a 
direction quite different from the deep-layer flow. In this case, the control analysis 
approximated the shallow-layer flow a bit more closely than the experimental analysis, and as 
a consequence the model control forecast without the winds was closer to the observed storm 
track. On a more positive note, Table 3 shows that thirteen of the eighteen 72 h forecasts were 
improvements over the control forecasts, with five of them being notable improvements 
(greater than 100 Km). Three of the notably improved forecasts occurred during Hurricane 
Andrew's interaction with an upper-level cyclonic circulation. An example is shown in Fig. 3, 
from the model initiation time of 1200 UT on 19 August 1992. At this time, the upper-level 
circulation was situated just to the north of Andrew and quite evident in the water vapor image 
sequences. This feature was well-described by the derived water vapor winds. During the next 
72 h, Fig. 3 shows that Andrew's actual track was only slightly affected by die upper-level low, 
shifting it northwest for a short time, before Andrew escaped its influence and turned more to 
the west. However, the VICBAR control forecast from this time recurved Andrew to the north 
and eventually to the northeast in response to the influence of the upper low on the deep layer 
mean steering flow. Although the turn back towards the west was not predicted, the 
experimental VICBAR forecast made with water vapor wind data responded with much less 
curvature, and an improved longer-range forecast. 
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b) Impact on regional numerical forecasts 


As part of a benchmark product demonstration, water vapor wind sets were produced twice 
a day at CIMSS during a six-week period in March-April 1994. Full disk wind sets (Fig. 2) were 
derived from both GOES-7 and Meteosat-3 around 00 and 12 UT each day. Processing was 
accomplished in a manner to simulate real time operations, and the winds were electronically 
delivered to NMC files in a time frame commensurate to meet operational deadlines. 

It is of interest to examine the impact of the water vapor winds on regional scale numerical 
analysis and prediction. To accomplish this, forecast impact tests were conducted employing 
the CIMSS regional assimilation system (CRAS), which is an adaptation of the Australian 
Bureau of Meteorology operational regional forecast system (Leslie et al. 1985). The model 
analysis domain covers much of North America and adjacent waters at 150 km horizontal 
resolution (Fig. 4), with 39 vertical levels. The prediction model operates with 20 vertical levels, 
and produces forecasts out to 72 h. Within the analysis package, a 3-dimensional variational 
wind-mass adjustment routine is activated in areas lacking in mass information (i.e. oceanic 
regions). This routine essentially attempts to exploit and spread the relative wealth of 
differential information (gradients implied by the satellite wind information) onto the mass field 
through gradient balance adjustments (Seaman et al. 1977). 

Experimental forecasts (EXP) which include the water vapor winds will be compared with 
control (CON) runs which initialized on operationally available data only. It should be noted that 
this exercise was not conducted in an assimilation mode; all analyses were independently 
produced with operational (NMC global model forecasts) boundary conditions only. An 
example of data impact on the CRAS analysis and forecast fields is illustrated in Figs. 4 and 5, 
respectively. Shown in Fig. 4 is a plot of the 275-325 Mb water vapor winds (large vectors) for 
00 UT on 22 March 1994. It should be noted that winds were intentionally produced only over 
marine areas. The vectors are overlain on a 300 Mb model wind analysis difference field (EXP 
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- CON), indicated by the smaller barbs (a barb length of one grid point interval equals 8 m/sec). 
In general, differences at this particular level are modest, but can be as great as 6 m/sec such 
as in the trough off of the east coast of the United States. 

These initial analysis differences translate into divergent 48 h predictions, illustrated as a 
300 Mb wind forecast difference field (EXP - CON) in Fig. 5. The 48 h EXP forecast 300 Mb 
height field is overlain to show the large-scale pattern. In order to highlight the differences, in 
this figure a barb length of one grid point interval equals 5 m/sec. Rg. 5 shows that while small 
differences cover most of the domain, more notable differences are found in several locations 
including the EXP forecast tendency to amplify the ridge over the central United States. 

In general, overall model impact results from CRAS forecasts of 52 cases from the 
March/April 1994 trial period reveal a modest but positive improvement in RMS error and bias 
statistics of upper-level height and wind fields. An example of the 400 Mb height forecast 
impact (72 h) in terms of RMS errors is given in Fig. 6. In this figure, a comparison of CON and 
EXP forecast RMS errors (versus verifying rawinsondes) is presented as a CON minus EXP 
value (meters) for each case. Positive model forecast impacts (reduced RMS errors relative to 
CON as indicated by cases above the zero line) produced by the addition of the water vapor 
winds to the analyses are indicated in a large majority of the 52 cases. These results are fairly 
representative for mandatory levels between 200 and 500 Mb, where the greatest impact 
would be expected due to the vertical distribution of the data. 

The impact of the water vapor winds on the regional model forecasts appears to be highly 
case dependent. While the absolute forecast error improvements are not large in a domain- 
wide RMS error sense, sub-regional forecasts may be more significantly affected. Qualitative 
examination of the cases run with the CRAS model indicates that the greatest positive impact 
occurs with active patterns in the eastern Pacific Ocean (i.e. upper-level circulations or jet 
streaks). These active periods correspond to some of the positive impact peaks in Rg. 6. Often 
these systems are not adequately defined by operationally available data (i.e. aircraft reports 
or cloud tracked winds). For an example, see Velden et al. (1994). 
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c. Impact on global numerical prediction 


As mentioned earlier, the water vapor wind sets produced by CIMSS for the trial period in 
March-April 1994 were transmitted to NMC for global model impact tests. NMC conducted the 
exercise using a T62 version of the global spectral model (similar to the operational T126 
AVN/MRF). The analysis component of this system is referred to as Spectral Statistical 
Interpolation (SSI), which differs only slightly from the previously employed optimum 
interpolation (01). Spectral coefficients used in the forecast model are analyzed directly using 
the same basic equations as the 01. However, this method is considered superior because no 
data selection is necessary, all observations are used simultaneously, and the dynamical 
constraint between the wind and mass fields is more realistic and applied globally (Derber et 
al., 1991). 

Forecasts produced from analyses initialized with the water vapor winds are compared to 
operational forecasts in order to assess model impact. Shown in Table 4 are comparisons of 
results from 72 h forecasts for a 15-day subset of the trial period. The forecast fields are 
verified against model analyses. Two latitude-dependent domains are considered in order to 
examine regional impacts. 

Considering first the near-global domain, the impacts are not dramatic, as might be 
expected since the wind coverage over the mid latitudes is limited. It should also be noted that 
the wind sets produced by CIMSS cover only a portion of the western hemisphere, while the 
verification statistics are global. Oespite these caveats, the impacts are slightly positive in all 
parameters examined. Nine of the fifteen forecasts were improved, with a slight reduction in 
the overall 200 Mb wind vector RMS error (results were similar at other upper-levels). 

Greater impact is seen in the tropical band, a region of sparse conventional data, and 
abundant water vapor wind coverage. While still modest, a more substantial gain in forecast 
accuracy is noted, especially considering the rotational component of the wind (psi). For this 
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particular parameter, twelve of the fifteen forecasts were improved. These results suggest the 
winds are contributing to a better definition of circulations or wave-like features in the ITCZ 
region. 

4. DISCUSSION 

There is little doubt that upper-tropospheric wind information derived from motions deduced 
in sequential geostationary satellite water vapor imagery can compliment the current 
operational observational network, especially over marine (data void) areas. High-density, full- 
disk wind sets from multiple satellites are currently attainable, and near-global coverage will be 
realized after successful deployment of the GMS-5 satellite in 1995. The observations are 
complimentary to upper-tropospheric cloud-drift winds by providing information in cloud-free 
areas, and statistical evaluations show they are comparable in quality. The global data sets 
can reliably be produced at 6 h intervals to meet current operational requirements, but could 
potentially be derived at hourly intervals. 

The limited NWP impact experiments presented here are encouraging. Three different 
NWP systems are examined, ranging from global, to regional, to hurricane track forecast 
prediction systems in order to evaluate the impact of the wind information on various 
atmospheric scales. Trials (forecasts) were conducted based on inclusion of the water vapor 
winds into the initial analyses, versus control runs without the winds. While the results are 
modest at best, positive forecast impact is realized with all three NWP systems. 

The treatment and incorporation of the wind observations by the respective analysis (or 
initialization) schemes is of course of primary interest in terms of realizing numerical model 
forecast impact. The limited magnitude of the impact results presented above could be partially 
explained by the analysis schemes, or the methodology employed to evaluate the impact For 
example, the global model for our trial was run by NMC at a reduced resolution which could 
limit data impact In regards to the regional model (CRAS), the evaluation employed 
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independent static analyses, while a pure assimilation mode may yield greater impact. The 
question is also raised whether the variational mass-wind blending procedure employed within 
the CRAS analysis scheme is optimal. Finally, the analysis module of the hurricane track 
forecast system (VICBAR) contains no vertical blending of information above 350mb (where 
most of the water vapor wind vectors reside). Instead, the information is only assigned to the 
nearest mandatory pressure level. The VICBAR initializes on a mass-weighted deep layer- 
mean (DLM) wind analysis created from mandatory level wind analyses between 850 Mb and 
200 Mb. Due to the mass weighting, the upper levels contribute a relatively low influence on 
the DLM field. This element of the analysis system is crucial to our applications, since the 
water vapor wind information is not being vertically spread to lower levels that have a greater 
influence on the DLM (which approximates the hurricane steering currents). The somewhat 
limited domain and bogus vortex procedure of the VICBAR may also place inherent limitations 
on the data impact. Despite the potential limiting factors outlined above, it is encouraging that 
positive forecast impact from the water vapor wind information was realized in each NWP 
experiment. Optimal methods to integrate these types of observations (i.e. 4D VAR?) will be 
investigated in future trials. 

Also encouraging are the prospects of water vapor winds from the recently launched 
GOES-8 satellite, which contains an imager that is characterized by improved resolution and 
signal-to-noise. Preliminary investigations are showing an improved water vapor wind product 
relative to GOES-7/Meteosat-3 (Velden et al. 1995). Further testing and algorithm refinement is 
underway, and NWP impact tests with NMC are planned in 1995. Upon a successful 
demonstration, we expect the GOES-8 water vapor winds to become part of the operational 
data base sometime in 1995. 

The GOES-8 satellite also carries a sounding instrument, which contains several water 
vapor channels which can sample the atmosphere at various levels. Potentially, imagery from 
these channels can be sequenced to provide motion (wind) observations at lower tropospheric 
levels. Combined with the upper-level observations from the imager, this could provide 
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important wind 'profile’ information. A candidate application would be hurricane steering 
currents, which are greatly influenced by the mid-tropospheric flow. This potential source of 
wind data is being investigated. 

Corresponding author address: Christopher Velden, Univ. of Wisconsin-SSEC, 1225 W. 
Dayton St., Madison, Wis., 53706. 


Acknowledgments: The author would like to express gratitude to Steve Nieman and Steve 
Wanzong of CIMSS, Robert Aune and Kit Hayden of NOAA/NESDIS, James Franklin and Sim 
Aberson of AOML/HRD, Steve Lord and Pete Kaplan of NMC, and Paul Menzel of 
NOAA/NESDIS. This research is being supported by NASA grant NAG8-974 and NOAA grant 
50 WC N E-30607 5 . 


References 

DeMaria, M., S.D. Aberson and K.V. Ooyama, 1992: A nested spectral model for hurricane 
track forecasting. Mon. Wea. Rev., 120 , 1628-1643. 

Derber, J.C., D. F. Parrish and S.J. Lord, 1991: The new global operational analysis system at 
the National Meteorological Center. Wea. and Forecasting, 6, 538-547. 

Franklin, J.L and M. DeMaria, 1992: The impact of Omega dropwindsonde observations on 
barotropic hurricane track forecasts. Mon. Wea. Rev., 120 , 381-391. 


14 


Hayden, C M. and R.J. Purser, 1988: Three-dimensional recursive filter objective analysis of 
meteorological fields. 8th Co nf. Numerical Weather Prediction, Baltimore, 185-190. 

Hayden, C.M. and T.R. Stewart, 1987: An update on cloud and water vapor tracers for 
providing wind estimates. 6th Symposium Met. Obs. and Instr., New Orleans, 70-75. 

Hayden, C.M. and C.S. Velden, 1991: Quality control and assimilation experiments with 
satellite-derived wind estimates. 9th Conf. Numerical Weather Prediction, Denver, 19-23. 

Holmlund, K., 1993: Operational water vapor wind vectors from Meteosat imagery. 2nd 
workshop on wind extraction from operational satellite data, Tokyo, 77-84. Available from 
EUMETSAT, Darmstadt, Germany. 

Laurent, H., 1993: Wind extraction from Meteosat water vapor channel image data. J. Appl. 
Meteor., 32, 1124-1133. 

Leslie, L.M., G.A. Mills, L.W. Logan, D.J. Gauntlett, G.A. Kelly, MJ. Manton, J.L. McGregor 
and J. Sardie, 1985: A high resolution primitive equations model for operations and research. 
Aust. Met. Mag., 33, 11-35. 

Merrill, R.T., W.P. Menzel, W. Baker, J. Lynch and E. Legg, 1991: A report on the recent 
demonstration of NOAA's upgraded capability to derive cloud motion satellite winds. Bull. 
Amer. Meteor. Soc., 72, 372-376. 

Nieman, S.J., W.P. Menzel, C.M. Hayden, S. Wanzong, and C.S. Velden, 1995: The NESDIS 
automated satellite derived winds system. Submitted to Bull. Amer. Meteor. Soc. 


15 


Seaman, R.S., R.L. Falconer and J. Brown, 1977: Application of a variational blending 
technique to numerical analysis in the Australian region. Aust. Met. Mag., 25, 3-24. 

Stewart, T.R., C.M. Hayden and W.L. Smith, 1985: A note on water-vapor wind tracking using 
VAS data on MclDAS. Bull. Amer. Meteor. Soc., 66, 1111-1115. 

Velden, C.S., C.M. Hayden, W.P. Menzel, J.L. Franklin and J.S. Lynch, 1992: The impact of 
satellite-derived winds on numerical hurricane track forecasting. Wea. and Forecasting, 7, 107- 
118. 


Velden, C.S., 1993: Investigation of water vapor winds from geostationary satellites. 2nd 
workshop on wind extraction from operational satellite data, Tokyo, 99-104. Available from 
EUMETSAT, Darmstadt, Germany. 

Velden, C.S., S. Nieman and S. Wanzong, 1994: Investigation of water vapor motion winds 
from geostationary satellites. 7th Conf. Satellite Meteor., Monterey, 360-363. 

Velden, C.S., 1995: Impact of geostationary satellite water vapor channel data on weather 
analysis and forecasting. 14th Conf. Wea. Analysis and Fore., Dallas, xxx-xxx. 

Weldon, R.J., and S.J. Holmes, 1991: Water vapor imagery: Interpretation and applications to 
weather analysis and forecasting. NOAA Tech. Report NESDIS 57, 5200 Auth Rd, Wash., 
D.C., 213pp. 


16 


List of Figures 


Figure 1 . Schematic of the objective quality control (QC) scheme employed with automated 
water vapor motion vectors (WVMV) production at CIMSS. The three-step process is 
designated by the numbers shown. 

Figure 2. Example of the horizontal coverage that can be attained from combined GOES-7 and 
Meteosat-3 water vapor wind data sets. 

Figure 3. Example of the impact of water vapor wind data (EXP) on a numerical track forecast 
(VICBAR) of Hurricane Andrew on 19 August 1992. The hurricane symbols represent the best 
track verification positions at 12, 24, 36, 48 and 72 hrs after initialization, and the CON 
(control) track is the operational forecast. 

Figure 4. Plot of 275-325 Mb water vapor winds (large barbs, full stick = 5 m/sec and flag = 25 
m/sec) over a 300 Mb CRAS model wind analysis difference (EXP minus CON) field (arrow 
length of one grid length interval = 8 m/sec) for 00 UT 22 March 1994. 

Figure 5. Difference field (EXP minus CON) from two CRAS model 48 h forecasts of 300 Mb 
wind (arrow length of one grid interval - 5 m/sec) initialized at 00 UT 22 March 1994. Also 
shown (contours in meters) is the attending 300 Mb height forecast. 

Figure 6. Comparison of RMS errors (verified against rawinsondes) for CRAS model 72 h 
forecasts of 400 Mb heights for 52 cases in March/April 1994, expressed in terms of 
RMS(CON) minus RMS(EXP). Positive impact of the water vapor wind data (EXP) is indicated 
in cases with points above the zero line. 
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Table 1. Impact of Meteosat-3 water vapor winds on VICBAR tropical cyclone track forecasts. 
The following verification statistics are valid for forecasts which included the water vapor 
winds: number of forecasts (N), effective number of independent forecasts (N*), and the mean 
forecast error (MFE) relative to best track verification. Also given are comparisons with the 
control forecasts: the mean forecast improvement (MFI), expressed in both kilometers and as a 
percent relative to the control forecast error, the standard deviation of the improvements (SOI), 
the number of improved forecasts (IF), the frequency of improved forecasts (FIF) expressed as 
a percent, and whether the forecast improvements are statistically significant at the 95% 
confidence level (SIG). 

Forecast 


Interval 

<h) 

N 

N* 

MFE 

(km) 

MFI 

(km) 

MFI 

% 

SOI 

(km) 

IF 

# 

FIF 

% 

SIG 

(y/n) 

24 

32 

26.6 

149.1 

3.3 

2.2 

23.8 

19 

60 

n 

48 

27 

22.4 

338.4 

3.9 

1.1 

48.2 

16 

60 

n 

72 

21 

17.4 

347.0 

23.6 

6.8 

91.6 

15 

71 

n 


Table 2. Same as Table 1 , except only for cases east of 70 W. 


Forecast 


Interval 

(h) 

N 

N* 

MFE 

(km) 

MFI 

(km) 

MFI 

% 

SOI 

(km) 

IF 

# 

FIF 

% 

SIG 

(y/n) 

24 

23 

19 

136.4 

9.0 

6.6 

20.2 

16 

70 

n 

48 

20 

17.8 

335.3 

10.2 

4.3 

60.1 

13 

65 

n 

72 

18 

16.2 

362.3 

30.4 

8.4 

98.6 

13 

72 

n 
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Table 3. Distribution of 72 h experimental forecast errors relative to control forecast 
errors for cases east of 70 W. 

Forecast differences (km) 

• ************** * * * 


200 


100 0 100 200 
BETTER WORSE 


Table 4. Impact of water vapor motion wind data 
on 72 h forecasts from the NMC global spectral 
model during a 15-day trial period in March, 1994. 

60 N - 60 S 





fsp 


opn 

exp 

exp/opn 

200 Mb vector rms (m/s) 

12.00 

11.86 

9/6 

200 Mb acpsi 

.809 

.815 

9/6 

200 Mb acv 

.741 

.742 

8/7 

20 N- 

20 S 



200 Mb vector rms 

10.51 

10.22 

9/6 

200 Mb acpsi 

.654 

.681 

12/3 

200 Mb acv 

.541 

.550 

8/7 


opn = operational runs 

exp = experimental runs with water vapor winds 
fsp = frequency of superior predictions 
acpsi = anomaly correlation of stream function 
acv = anomaly correlation of v-component 
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